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Preface 


In  order  to  more  accurately  predict  nuclear  near 
and  far  field  fallout  effects,  an  accurate  portrayal  of  the 
stabilized  nuclear  cloud  particle  size  distribution  is  a 
necessity.  However,  there  is  a  continuing  debate  over  the 
nature  of  this  distribution.  During  the  late  1950s  and 
1960s  much  research  was  done  using  ground  fallout  samples 
to  construct  the  nuclear  cloud  particle  size  distribution. 
Unfortunately,  ground  fallout  samples  cannot  adequately 
describe  the  particle  size  distribution  in  the  stabilized 
cloud  and  airborne  sampling  of  our  early  events  was  scarce 
and  poorly  documented.  Consequently,  in  most  cases,  the 
nuclear  cloud  particle  size  distribution  currently  being 
used  by  defense  planners  to  predict  fallout  patterns  of 
strategic  importance  is  based  on  Marcel  Nathans'  1970  pub¬ 
lication  of  his  work  with  nuclear  cloud  samples.  Yet, 
aside  from  Nathans'  work,  little  research  has  been  done  to 
derive  a  more  accurate  nuclear  cloud  particle  size  distribu 
tion  from  nuclear  cloud  samples. 

The  purpose  of  this  independent  study  was  to  try 
a  different  approach  to  reconstructing  the  stabilized 
nuclear  cloud  particle  size  distribution  from  reduced  air¬ 
borne  filter  sample  data.  This  method  requires  only  one 
piece  of  reduced  filter  sample  data,  the  total  mass  or  the 


total  activity  on  the  filter.  It  also  uses  a  state-of- 
the-art  stabilized  cloud  model  and  state-of-the-art 
particle  fall  mechanics. 

Overall,  this  research  was  quite  enjoyable.  I 
would  like  to  thank  Dr.  Charles  J.  Bridgman  for  his  sup¬ 
port  and  patience  with  me.  Also,  I  am  deeply  indebted  to 
my  wife,  4PM/  f°r  her  Pafcience  an^  understanding  during 
the  many  hours  needed  to  complete  this  work. 
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Abstract 


This  study  developed  a  numerical  method  of  unfold¬ 
ing  the  particle  size  distribution  of  the  stabilized 
nuclear  cloud  from  reduced  airborne  filter  sample  data^  ~ / 
A  stabilized  nuclear  cloud  is  modeled  using  a  trial  par¬ 
ticle  size  distribution  that  is  positioned  in  the  atmos¬ 
phere  by  empirical  relationships . developed  by  Hopkins  and 
Connors ^  Davies-McDonald  fall  mechanics  are  used  to  model 
the  falling  particles  in  the  ‘iuoieat,  cloud.  The  amount  of 
mass  at  each  sample  altitude,  at  each  sample  time  is  calcu¬ 
lated  ‘from  the  cloud  model  and  compared  to  the  amount  of 
mass  found  in  4h^  actual  cloud  samples.  When  the  calcu¬ 
lated  masses  equal  the  actual  masses,  the  particle  distri¬ 
bution  used  to  construct  the  stabilized  cloud  is  the 
correct  one.  A  computer  code  for  this  numerical  analysis 
is  also  presented. 

The  computer  code  is  tested  using  hypothetical 
filter  sample  data  constructed  from  a  known  particle  size 
distribution.  Additionally ,>  an  input  parameter  sensitivity 

V<' 

analysis  is  conducted. 

Actual  nuclear  cloud  sample  data  from  the  Redwing 
series,  shot  ZUNI  is  analyzed  using  this  numerical  method 
of  airborne  nuclear  cloud  sample  analysis.  The  outcome  of 
the  ZUNI  sample  analysis  is  somewhat  inconclusive  in  that 


IX 


it  does  not  pinpoint  a  distribution.  However,  based  on 

"  > 

the  results  of  the  model  sensitivity  analysis,  the  ZUNI 
sample  analysis  indicates  that  the  particle  size  distribu 
tion  of  the  stabilized  ZUNI  cloud  may  be  lognormal  with  a 
log-slope  that  i§> between  2.9  and  3. 9* but  is  definitely 
not  less  than  2.7  nor  greater  than  5.0. 

Finally,  the  theory  for  an  alternative  method  of 
airborne  sample  analysis  is  presented.  This  method  uses 
the  relative  number  of  particles  of  each  size  found  in  an 
airborne  sample  to  unfold  the  stabilized  nuclear  cloud 
particle  size  distribution. 


A  NUMERICAL  METHOD  FOR  UNFOLDING  THE  STABILIZED 


NUCLEAR  CLOUD  PARTICLE  DISTRIBUTION 


I.  Introduction 


Background 

The  ability  to  accurately  model  fallout  patterns 
from  nuclear  bursts  is  largely  dependent  upon  knowledge  of 
the  particle  size  distribution  in  the  stabilized  nuclear 
cloud.  Bigelow’s  sensitivity  analysis  of  lognormal  par¬ 
ticle  size  distributions  versus  fallout  pattern  predictions 
(3:V-15)  graphically  illustrates  this  fact.  Moreover,  he 
demonstrated  that  selection  of  the  proper  distribution 
standard  deviation  is  very  important  because  relatively 
small  changes  in  this  value  produce  changes  of  an  order  of 
magnitude  or  more  in  the  size-activity  distribution  median. 
These  size-activity  median  changes  greatly  affect  the  fall¬ 
out  prediction  outcome  (3:V-13).  Consequently,  the  impor¬ 
tance  of  selecting  the  proper  nuclear  cloud  particle  size 
distribution  for  fallout  pattern  modeling  cannot  be  empha¬ 
sized  strongly  enough. 

Many  proposed  nuclear  cloud  particle  size  distribu¬ 
tions  can  be  found  in  the  open  literature.  For  the  most 
part,  these  distributions  were  determined  by  evaluation  of 
fallout  samples  from  the  many  nuclear  test  shots  in  Nevada 


and  the  Pacific  ocean.  Conners  lists  twelve  such  distribu¬ 
tions  in  his  study  (7:11).  Two  of  the  most  popular  ones 
are  the  Defense  Land  Fallout  Interpretive  Code  (DELFIC) 
default  distribution  and  the  so-called  TTAPS  distribution. 
These  two  lognormal  particle  size  distributions  are  radi¬ 
cally  different,  yet  they  were  both  derived  from  fallout 
sample  data.  Therefore,  it  is  quite  clear  that  there  is 
still  a  great  potential  for  debate  about  the  nature  of  the 
stabilized  nuclear  cloud  particle  size  distribution. 

Another  area  of  growing  concern  is  over  the  poten¬ 
tial  differences  in  the  particle  size  distribution  calcu¬ 
lated  from  analysis  of  "fallen"  fallout  samples  and  the 
actual  "falling"  particle  size  distribution  in  the  nuclear 
cloud.  Hopkins'  analysis  of  the  Mount  St.  Helens  ash  cloud 
determined  that  the  particle  size  distribution  calculated 
from  "fallen"  samples,  collected  on  the  ground,  could  not 
be  used  to  adequately  predict  the  observed  ash  fallout  pat¬ 
tern  (16:83-91).  It  is  quite  possible  that  this  phenomenon 
is  true  of  nuclear  clouds  as  well. 

Overall,  the  main  issue  here  is  that  there  is  uncer¬ 
tainty  surrounding  every  potential  particle  distribution 
used  to  describe  the  stabilized  nuclear  cloud.  Since  the 
selection  of  a  proper  particle  size  distribution  is  key  to 
accurate  fallout  prediction  it  is  imperative  that  every 
effort  be  made  to  accurately  define  a  valid  distribution. 
Therefore  it  is  evident  that  analysis  of  "fallen"  particle 


samples,  taken  on  the  ground,  must  be  used  in  conjunction 
with  analysis  of  "falling"  samples,  taken  in  the  cloud,  in 
order  to  better  qualify  the  actual  particle  size  distribu¬ 
tion  in  the  nuclear  cloud.  This  study  is  primarily  con¬ 
cerned  with  the  analysis  of  "falling"  cloud  samples. 

Problem 

Much  effort  has  been  made  in  analyzing  the  "fallen" 
fallout  samples.  However,  aside  from  Nathans'  work  (Ref 
19) ,  little  has  been  done  to  unfold  the  nuclear  cloud  par¬ 
ticle  size  distribution  from  airborne  "falling"  samples. 
This  study  attempts  to  unfold  the  particle  size  distribu¬ 
tion  of  the  stabilized  nuclear  cloud  through  numerical 
analysis  of  airborne  filter  samples  taken  from  nuclear 
clouds  during  atmospheric  testing. 

Scope 

This  study  is  limited  in  scope  to  examination  of 
samples  taken  from  nuclear  clouds  by  aircraft.  The  model 
used  here  assumes  a  horizontal  flight  path  through  the 
center  of  the  nuclear  cloud  at  different  altitudes  and  dif¬ 
ferent  times  following  cloud  stabilization.  The  mass  of 
material  collected  on  the  filter  samples  is  then  compared 
to  what  should  have  been  on  the  filter  samples  given  a 
hypothetical  nuclear  cloud  particle  size  distribution. 

The  stabilized  nuclear  cloud  is  modeled  as  being 


gaussian  in  the  "x-y"  horizontal  plane  as  well  as  gaussian 
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in  the  vertical  or  "z"  direction.  The  particles  are 
gravity  sorted  by  size.  In  other  words,  each  size  group  is 
lofted  to  its  initial  stabilized  altitude  and  spatially 
positioned  using  Hopkins'  and  Conners'  empirical  relation¬ 
ships.  No  wind  effects  are  considered.  Only  the  forces  of 
gravity  and  atmospheric  viscosity  act  on  the  particles. 
Davies-McDonald  fall  mechanics  are  used  to  model  particles 
falling  through  a  non-homogeneous  atmosphere. 

Only  surface  bursts  or  near  surface  bursts  are  con¬ 
sidered  in  this  study.  Surface  bursts  loft  the  most 
material  and  their  nuclear  clouds  are  believed  to  be  better 
understood  than  clouds  from  near  surface  or  air  bursts. 

Finally,  the  lognormal  distribution  is  the  only 
distribution  considered  as  a  potential  nuclear  cloud  par¬ 
ticle  size  distribution  in  this  study.  It  is  hypothesized 
that  a  reasonable  fit  can  be  achieved  with  a  lognormal  dis¬ 
tribution  (4:209).  Moreover,  numerical  analysis  of  either 
the  mass  or  activity  distributions  in  the  nuclear  cloud  is 
greatly  simplified  by  the  property  of  lognormal  moments  for 
a  lognormal  distribution. 

Assumptions 

Minor  assumptions  concerning  specific  procedures 
used  in  this  study  will  be  presented  at  appropriate  places 
in  the  body  of  this  report.  However,  the  following  are 
some  general  assumptions  that  apply  throughout  the  work. 
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1.  The  stabilized  nuclear  cloud  can  be  represented 
by  a  gaussian  distribution  of  all  particle  sizes  in  the 
"x-y"  horizontal  plane.  In  the  vertical  or  "z"  direction 
the  nuclear  cloud  can  be  modeled  by  using  a  finite  number 
of  particle  size  groups.  Each  group  is  represented  by  a 
mean  radius  and  normally  distributed  in  the  vertical  direc¬ 
tion  . 

2.  There  is  homogeneous  mixing  of  particles  hori¬ 
zontally  within  the  cloud. 

3.  The  particle  size  distribution  in  the  stabil¬ 
ized  nuclear  cloud  can  be  modeled  with  a  lognormal  distri¬ 
bution  . 

4.  Gravity  and  atmospheric  viscosity  are  the  only 
forces  acting  on  the  particles  in  the  stabilized  nuclear 
cloud.  No  wind  effects  are  considered. 

5.  All  samples  were  taken  under  the  same  condi¬ 
tions;  the  same  sampling  apparatus  was  used  for  each  sample 
Additionally,  the  aircraft  taking  the  samples  was  flown 
through  the  horizontal  geometric  center  of  the  cloud  at 
each  respective  sample  altitude. 

6.  All  fallout  particles  are  assumed  to  be  spheri¬ 
cal.  This  assumption  has  little  effect  on  fallout  predic¬ 
tion  (21:32) . 

7.  The  fallout  density  is  assumed  to  be  constant, 
2600  kilograms  per  cubic  meter,  independent  of  particle 


roach  and  Sequence  of  Presentation 

The  computer  model  used  for  airborne  filter  sample 
analysis  is  developed  in  Chapter  II  in  the  following 
sequence.  The  stabilized  nuclear  cloud  and  falling  cloud 
models  are  presented  first.  Next,  the  particle  distribu¬ 
tion  calculative  method  is  presented.  Then,  nuclear  cloud 
sampling  theory  and  integration  of  sampling  into  the  fall¬ 
ing  cloud  model  are  described.  Finally,  the  method  for 
determining  the  particle  size  distribution  through  falling 
particle  sedimentation  and  numerical  analysis  of  reduced 
filter  sample  data  is  presented  in  algorithm  form. 

In  Chapter  III,  the  model  is  validated  by  using 
calculated  hypothetical  filter  samples  from  a  known  par¬ 
ticle  size  distribution  as  input  to  the  computer  code  pre¬ 
sented  in  Chapter  II.  This  code  unfolds  the  stabilized 
nuclear  cloud  particle  size  distribution  from  the  input 
filter  sample  data.  Additionally,  model  sensitivity  to 
input  variations  is  presented  and  discussed  in  the  final 
section  of  this  chapter. 

In  Chapter  IV,  a  study  of  actual  airborne  sample 
data  from  the  Redwing  series,  shot  ZUNI,  is  presented  with 
results . 

Finally,  an  alternative  method  of  reconstructing  the 
stabilized  nuclear  cloud  particle  size  distribution  from  one 
or  two  airborne  filter  samples  is  presented  in  Chapter  V. 

Overall  results  and  conclusions  are  found  in  Chap- 


II.  Model  Development  and  Theory 


Background 

Nuclear  cloud  particle  formation  is  a  complicated 
process  that  occurs  during  the  time  that  the  nuclear  cloud 
is  cooling.  Initially,  the  nuclear  detonation  releases 
massive  quantities  of  energy  in  a  very  short  period  of  time. 
The  resultant  high  temperature  from  x-ray  deposition  in  the 
atmosphere  creates  a  fireball  that  literally  vaporizes 
everything  within  its  boundaries.  These  vaporized  materials 
include  unfissioned  weapons  material,  fission  fragments, 
weapons  case  material,  and  any  soil  that  was  consumed  by 
the  fireball  in  the  case  of  a  surface  or  near-surface  burst. 
As  the  fireball  expands  and  rises,  it  cools  and  the  highly 
refractory  vaporized  material  begins  to  condense  and  form 
particles.  Additional  material  may  be  transported  up  into 
the  hot  rising  cloud  by  updrafts  created  by  the  rapidly 
rising  hot  cloud.  Generally,  much  of  this  material  is  not 
completely  vaporized.  Once  in  the  cloud,  these  particles 
grow  in  size  as  the  more  volatile  elements,  still  in  a 
vapor  state,  plate-out  on  their  surfaces.  The  volatile 
material  continues  to  plate-out  on  the  surfaces  of  all  par¬ 
ticles  in  the  cloud  until  the  cloud  temperature  achieves 
an  equilibrium  with  ambient  temperatures. 


Stabilization  occurs  when  the  hot  nuclear  cloud  has 
cooled  to  the  point  where  it  stops  rising.  Regardless  of 
yield,  stabilization  usually  occurs  within  5  to  7  minutes 
after  the  detonation.  However,  the  greater  the  weapon 
yield,  the  higher  its  cloud  will  rise  prior  to  stabiliza¬ 
tion.  For  example,  the  cloud  from  a  10  kiloton  weapon, 
surface  burst,  will  stabilize  at  approximately  5000  meters, 
whereas  the  cloud  from  a  1  megaton  surface  burst  will 
stabilize  at  approximately  13000  meters  (14:431).  Since 
the  scope  of  this  study  is  limited  to  times  following  cloud 
stabilization,  and  is  primarily  concerned  with  the  vertical 
distribution  of  particles  in  the  cloud,  the  results  of  a 
cloud  rise  model  to  be  used  in  this  study  are  described  in 
the  next  section. 

Stabilized  Cloud  Model 

The  altitude  to  which  a  particle  will  rise  depends 
on  the  weapon  yield  and  the  particle  mass.  For  example, 
the  smaller,  less  massive  particles  are  lofted  higher  by 
the  hot  rising  cloud  than  the  larger  more  massive  particles. 
Additionally,  because  of  the  assumption  that  the  nuclear 
cloud  is  gaussian  in  the  vertical  direction,  all  particles 
of  a  given  size  are  assumed  to  be  normally  distributed  in 
the  vertical  direction.  Hence,  the  gravitational  sorting 
and  particle  vertical  distribution  assumptions  lead  to  a 
cloud  rise  algorithm  for  modeling  the  stabilized  cloud. 

This  algorithm  uses  empirical  relationships,  presented  in 


Appendix  A,  to  calculate  the  lofted  altitude  of  the  center 
of  a  vertical  gaussian  distribution  and  standard  deviation 
for  each  different  particle  size  being  considered.  In 
reality,  there  are  a  great  number  of  particle  sizes,  there¬ 
fore,  an  equally  large  number  of  vertical  gaussian  distribu 
tions  to  be  lofted  and  characterized.  However,  this  study 
initially  limits  the  number  of  particle  size  groups  used  to 
at  least  the  number  of  cloud  samples  being  evaluated  and 
at  most  eight. 

An  illustrative  conception  of  the  vertical  distribu 
tion  in  the  stabilized  cloud  is  found  in  Figure  1.  A  gaus¬ 
sian  for  each  of  three  particle  size  distributions  is  shown 
The  highest  gaussian  represents  the  smallest  particles. 

Particle  Fall  Mechanics 

The  force  of  gravity  eventually  overcomes  the 
forces  that  are  causing  the  particles  to  rise  and  they 
begin  to  fall  from  the  cloud.  The  spherical  particle 
assumption  simplifies  the  problem  and  allows  a  particle's 
fall  velocity  to  be  calculated  based  on  its  radius  and 
altitude.  This  process  would  be  further  simplified  if  all 
particle  radii  were  less  than  about  10  micrometers  because 
Stokes  law  applies  to  particles  with  radii  less  than  10 
micrometers  falling  through  the  atmosphere.  However,  for 
particles  larger  than  10  micrometers,  aerodynamic  drag 
must  be  considered.  Therefore,  the  following  balance  of 
forces  equation  applies  (4:212): 
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Fig.  1.  Stabilized  Cloud  Particle  Size 
Spatial  Representation 
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where 

Pa  =  air  density  at  the  particle's  altitude  in 
kilograms/cubic  meter 

v  =  particle  velocity  in  meters/second 

Cd  =  drag  coefficient 

r  =  particle  radius  in  meters 

p^  =  particle  density  in  kilograms/cubic  meter 

g  =  gravitational  constant  in  meters/second 
squared 


mmm 


Since  the  drag  coefficient  is  also  a  function  of 
particle  velocity.  Equation  (1)  cannot  simply  be  solved 
for  the  particle  velocity.  However,  the  Reynolds  number 


for  falling  spheres  is  (18:464): 

2vpar 


(2) 


where 


v  =  particle  velocity  in  meters/ second 

p  =  air  density  at  the  particle's  altitude  in 
a  kilograms/cubic  meter 

r  =  particle  radius  in  meters 

n  =  dynamic  viscosity  of  the  air  at  the  particle's 
altitude  in  kilograms/meter-second 


Equation  (2)  can  be  solved  for  velocity  and  sub¬ 
stituted  into  Equation  (1) .  With  some  algebraic  manipula¬ 
tion,  this  yields: 


32p  p.  gr" 

r  cd  -  -.V 


3n 


(3) 
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Davies  (8:  259-270)  related  the  quantity  R  Cd  to  the 
Reynolds  number  by  the  following  two  empirical  relation¬ 
ships  : 


R 


-  2.3363  xlO"4  (R2Cd)2  +  2.0154  x  10 
x  (R2Cd)3  -  6.9105 x  10-9 (R2Cd)4 


LOG1Q(R)  =  -1.29536  +  0 .986  LOG1q (R2Cd)  -  0.046677 

x  [LOG1q (R2Cd) ] 2 + 0.0011235  [LOG1Q (R2Cd> ] 3 

for  3  <  R  <  10000  (5) 

The  Reynolds  number,  calculated  from  these  empirical  rela¬ 
tionships,  can  be  used  in  Equation  (2)  to  solve  for  the 
particle  velocity. 

Finally,  a  slip-drag  correction  factor  (22:6) 

SD  =  1  +  1.165  x  10  If rp  (6) 

a 

is  used  to  correct  the  falling  velocities  of  the  particles 
at  high  altitudes  for  their  reduced  interactions  with  air 
molecules  (4:212).  The  variables  r  and  p  are  the  same  as 

Cl 

defined  in  Equation  (2) . 

Now  that  the  fall  velocity  can  be  calculated  for 
any  particle  based  on  its  radius  and  altitude,  the  distance 
traveled  by  a  particle,  DZ ,  is  simply 

DZ  =  v  *  At  *  SD  (7) 

where  At  is  the  time  increment  in  seconds  for  which  the 
particle  falls. 

The  following  algorithm,  based  on  the  above  discus¬ 
sion,  is  used  in  this  study  to  model  particle  gravitational 


1.  U.S.  Standard  Atmosphere  is  used  to  calculate 


the  air  density  and  dynamic  viscosity  for  each  particle 
altitude . 

2 

2.  R  is  calculated  using  Equation  (3). 

2 

3.  With  R  the  Reynolds  number  is  calculated 
using  Equation  (4)  or  Equation  (5). 

4.  The  Reynolds  number  is  used  in  Equation  (2)  to 
determine  the  particle  velocity. 

5.  The  particle  velocity  is  corrected  for  slip- 
drag  by  Equation  (6) . 

6.  The  distance  fallen  is  calculated  using  Equa¬ 
tion  (7)  given  an  increment  of  fall  time.  This  distance 
is  subtracted  from  the  particle’s  altitude  yielding  a  new 
particle  altitude. 

7.  If  the  particle  has  fallen  for  the  desired 
time,  the  sequence  is  stopped;  if  not,  the  sequence  is 
repeated  for  additional  time  increments  until  the  total 
desired  fall  time  is  accomplished. 

The  optimum  time  increment  must  be  selected  so  that 
computer  time  is  conserved,  yet  sufficient  accuracy  is 
achieved.  Obviously,  selection  of  too  long  of  a  time  incre¬ 
ment  results  in  inaccuracies  because  of  variations  in 
atmospheric  properties  with  decreasing  altitude.  However, 
selection  of  too  short  of  an  increment  results  in  a  super¬ 
fluous  number  of  iterations  with  no  notable  increase  in 
accuracy.  Conners  (7:120)  suggested  that  a  time  increment 

13 


:*;»'*.* ‘I*,  ii.vw'v: 


be  selected  so  that  the  largest  particle  being  considered 
not  fall  more  than  1400  meters  during  that  time  increment. 
His  empirical  testing  determined  that  this  technique  pro¬ 
duced  the  largest  time  increment  that  did  not  significantly 
affect  the  calculated  distance  for  the  fall  of  particles. 
Therefore,  in  this  particle  fall  model,  the  time-of-fall 
increment  will  be  selected  based  on  the  time  it  takes  for 
the  largest  particle  size  being  considered  to  fall  1400 
meters . 

Lofted  Particle  Number 
Size  Distribution 

Given  that  all  the  previous  assumptions  are  plaus¬ 
ible,  one  final  assumption  is  needed  to  complete  the  model. 
The  lofted  particle  number-size  distribution  is  assumed 
to  be  lognormal  (4:210): 

.  2 

M  /In  r  -  a  \ 

N(r> - ^ -  .-*■( - S - ,8, 

/Tn  g  r 


where 


N(r)  =  the  number  of  particles  with  radius  r 

N  =  the  total  number  of  particles 

3  =  the  distribution  log-slope 

r  =  particle  radius  in  micrometers 

a  =  the  natural  log  of  the  distribution  median 
radius,  r 
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The  lognormal  distribution  is  a  universal  assumption  in 
approximating  nuclear  cloud  particle  size  distribution 
(12:7)  . 

Of  course,  one  very  useful  property  of  the  log¬ 
normal  distribution  is  that  the  moments  of  the  distributions 
are  also  lognormal  distributions  (2:12).  For  example, 
the  third  moment  of  the  lognormal  number-size  particle  dis¬ 
tribution  in  Equation  (8)  is  also  a  lognormal  distribution. 
This  is  the  mass-size  distribution: 

Mt  /ln  r  ~  a3\2 

M(r)  =  -  e  0  /  (9) 

/2ir  8  r 

where 

M(r)  =  the  total  mass  of  all  the  particles  with 
radius  r 

=  the  total  mass  lofted  by  the  weapon 

a.  =  a  +  3  *  (3^ 

3  o 

All  other  parameters  are  the  same  as  defined  above  follow¬ 
ing  Equation  (9) . 

When  integrated,  Equation  (9)  yields  the  total  mass 
aloft.  By  using  a  function  that  approximates  the  cumula¬ 
tive  lognormal  distribution,  the  mass  moment  can  be 
divided  into  any  number  of  equal  mass  groups  (1:932).  Then 
the  median  radius  for  each  equal  mass  group  can  be  used  to 
represent  the  entire  group  in  a  particle  fall  mechanics 
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Figure  2  also  graphically  illustrates  how  cloud 
samples  taken  at  different  altitudes  and  different  times 
intersect  the  various  particle  vertical  gaussian  distribu¬ 
tions  as  they  fall  with  time.  This  concept  is  the  key  to 
the  theory  behind  this  study's  analysis  of  airborne  samples. 
Simply  stated,  the  total  mass  of  material  in  an  airborne 
sample  is  equal  to  the  sum  of  the  contributions  of  mass 
from  each  equal  mass  group's  vertical  distribution  collected 
by  the  sampling  system  as  it  passed  through  the  nuclear 
cloud.  This  can  be  mathematically  stated: 

n 

s  =  yv  do) 

i=l 


where 

S  =  sample  total  mass 

n  =  number  of  equal  mass  groups  considered  (also 
equal  to  the  number  of  airborne  samples  being 
analyzed) 

G^  =  total  mass  contributed  by  each  equal  mass  group 


The  total  mass  contributed  by  each  group,  G^,  is  further 
defined  as: 


G.  =  A  *  F.  *  M, 
i  it 


(11) 


whei  j 


A  =  fraction  of  total  mass  aloft  contained  in 
group  i 


17 


F .  =  gaussian  fraction  of  group  i  at  the  sample 
1  altitude  and  time 

=  total  mass  lofted  by  the  weapon 

Note  that  if  the  correct  particle  mass  distribution  is  used, 
then  the  "A"  values  are  all  equal  since  equal  mass  groups 
were  assumed  initially. 

If  a  number  of  samples  are  taken  at  different  alti¬ 
tudes  and  at  different  times,  the  resulting  set  of  equa¬ 
tions  can  be  written: 


A-.F..M,  +  A_F.  +  ...  A  F.  ML  =  S. 
1  11  t  2  12  t  n  In  t  1 


AlF21Mt  +  A2F22Mt  +  ***  AnF2nMt  S2 


AlFnlMt  +  A2Fn2Mt  +  *  * •  AnFnnMt  =  Sn 


where,  in  the  case  of  these  equations, 

n  =  the  number  of  airborne  samples  being 
analyzed 

S..  .  .  =  samples,  each  taken  at  a  different 

altitude  and  time 

A,.  .  >  =  the  fraction  of  the  mass  aloft  con- 

(1  to  n)  ,  .  ,  .  , 

tamed  in  each  group 


F..  .  „  ,  .  =  the  gaussian  fraction  of  each  group 

at  each  sample  altitude,  at  each 
sample' time 

=  total  mass  lofted  by  the  weapon 


The  total  mass  on  "n"  airborne  samples,  or  "S" 


% 


values,  are  known.  By  using  the  stabilized  cloud  model 
described  above,  and  allowing  "n"  equal  mass  groups  to  fall 
in  accordance  with  the  fall  mechanics  described  above,  the 
^  gaussian  fraction,  or  "F"  value  of  each  of  the  equal  mass 

groups  can  be  calculated  for  each  sample  altitude  and  time. 
The  fractions  of  mass  aloft  contained  in  "n"  equal  mass 
^  groups,  or  "A"  values,  are  the  unknowns.  Consequently, 

this  type  of  analysis  yields  "n"  simultaneous  equations 
with  "n"  unknowns  that  can  be  expressed  in  a  simple  matrix 
&  equation: 


* 


m 


Mt[F]  X  [A]  =  [S]  (15) 

where 

F  -  n  x  n  matrix  (calculated) 

A  =  n  x  1  matrix  (unknowns) 

S  =  n  x  1  matrix  (sample  values) 

If  this  equation  is  multiplied  by  the  inverse  of  the  "F" 
matrix,  it  yields  a  solution  for  the  "A"  vector.  If  the 
solution  yields  equal  "A"  values,  then  the  fraction  of  the 
total  mass  lofted,  in  each  group,  is  the  same.  Since  the 
assumed  particle  distribution  was  initially  divided  into 
equal  mass  groups,  and  the  calculated  fraction  of  mass  per 
group  is  the  same,  then  the  assumed  particle  distribution 
is  the  correct  one. 
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The  solution  of  Equation  (15)  is  the  final  major 
step  in  determining  the  particle  size  distribution  for  the 
stabilized  nuclear  cloud.  From  now  on,  this  method  will  be 
referred  to  as  the  integral  method,  since  it  attempts  to 
unfold  the  stabilized  cloud  particle  distribution  by 
summing  the  contribution  from  falling  equal  mass  groups  to 
determine  the  total  sample  mass.  The  next  section  outlines 
the  code  used  to  do  the  distribution  search  thereby  inte¬ 
grating  all  the  previously  discussed  theory  into  a  distri¬ 
bution  search  algorithm. 

Airborne  Sample  Analysis  Code 
for  the  Integral  Method 

A  computer  program  was  developed  to  numerically 
analyze  airborne  cloud  sample  data  and  determine  a  par¬ 
ticle  distribution  for  the  stabilized  cloud.  This  program 
was  written  in  FORTRAN-5  (FORTRAN-77)  and  run  on  a  CYBER 
170-845.  The  program  is  called  SEARCH 4 .  Basically, 

SEARCH4  reads  a  log-slope  and  initial  median  radius  from  a 
data  file  called  SEARCH.  Then,  it  varies  the  median  radius 
through  a  range  of  values,  creating  a  number  of  different 
lognormal  distributions.  Each  distribution  is  checked  by 
the  algorithm  listed  below.  The  log-slopes  in  the  data 
file  SEARCH  vary  from  1.5  to  5.0.  The  median  radii  range 
of  values  is  selected  so  that  the  equal  mass  group  mean 
radii  calculated  from  each  lognormal  distribution  make 
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physical  sense.  A  copy  of  the  code,  example  input  files, 
and  glossary  of  program  variables  are  included  in  Appen¬ 
dix  B. 

The  algorithm  contained  in  the  program  is  as  fol¬ 
lows  : 

1.  Airborne  sample  data,  sample  total  mass,  time 
of  collection,  and  collection  altitude,  are  read  from  a 
file  "FILTER." 

2.  The  first  set  of  lognormal  distribution  param¬ 
eters  is  read  from  a  file  "SEARCH." 

3.  A  lognormal  distribution,  defined  by  a  log- 
slope,  g,  and  a  median  radius,  rQ,  is  divided  into  "n" 
equal  mass  groups  and  a  median  radius  for  each  of  the 
groups  is  calculated. 

4.  The  largest  median  radius  is  used  to  calculate 
the  fall  time  increment. 

5.  Each  group  represented  by  a  vertical  gaussian 
with  that  group's  median  radius  is  empirically  lofted  to 
its  stabilized  altitude  and  allowed  to  fall  until  the  time 
of  the  first  sample.  The  gaussian  is  evaluated  at  the 
first  sample  altitude  to  determine  the  fraction  of  total 


mass  in  this  group  being  contributed  to  the  total  sample 
mass.  The  mass  contributed  by  this  group  is  stored  in  an 
array,  . 
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6.  The  procedure  in  (5)  is  repeated  for  each 
group,  once  for  each  sample.  The  resultant  array,  F^ , 
is  the  "F"  matrix,  discussed  in  the  previous  section. 

7.  The  "F"  matrix  is  inverted  and  multipled  times 
the  "S"  vector  to  solve  Equation  (15)  for  the  fraction  of 
mass  in  each  group,  or  "A"  values. 

8.  The  "A"  values  are  compared  to  each  other  by 
least  squares  standard  deviation  and  only  those  with  a  low 
standard  deviation  are  written  to  the  output  file  "ANSWER" 
along  with  the  lognormal  distribution  parameters  for  that 
iteration. 

9 .  A  new  set  of  lognormal  distribution  parameters 
are  read  from  the  input  file  "SEARCH”  and  the  entire  pro¬ 
cess  is  repeated. 

In  summary,  here  is  an  example  that  will  demon¬ 
strate  the  integral  method  theorized  throughout  Chapter  II 
and  described  in  the  algorithm  in  this  section:  If  3  air¬ 
borne  sample  masses  are  to  be  analyzed,  then  a  trial  log¬ 
normal  mass-size  distribution,  characterized  by  a  log-slope 
and  a  median  radius,  is  divided  into  3  equal  mass  groups. 

A  median  radius  for  each  of  the  equal  mass  groups  is  cal¬ 
culated  and  used  to  represent  all  the  mass  in  its  group. 
Then,  each  of  the  equal  mass  groups,  with  median  radius, 
represented  by  a  vertical  gaussian,  is  allowed  to  fall 
from  its  stabilized  altitude  until  each  sample  time.  The 
vertical  gaussian  is  evaluated  at  the  sample  altitude  to 


mmSm&mm 

mmmm 

& 

"U wrr. 


'wimugnwiwnw^iTii Tirw 


determine  what  fraction  of  the  mass  represented  by  that 
group  is  contributed  to  the  total  mass  in  the  sample.  Now, 
all  chs  elements  for  Equation  (15)  are  present  and  the 
matrix  equation  can  be  solved  for  the  fraction  of  mass  in 
each  group,  or  the  "A"  values.  If  the  "A"  values  are  all 
equal,  then  the  trial  mass-size  distribution  is  the  correct 
one. 
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III.  Integral  Model  Validation 


Background 

Prior  to  analyzing  airborne  filter  sample  data  from 
an  actual  nuclear  shot,  the  integral  model  presented  in 
Chapter  II  was  validated  using  hypothetical  sample  data 
created  with  a  known  particle  distribution.  The  purpose  of 
this  validation  was  twofold.  First,  the  obvious  reason  for 
controlled  testing  of  the  model  was  to  insure  that  it  would 
function  as  predicted.  Second,  and  perhaps  of  equal  impor¬ 
tance,  was  the  need  to  conduct  a  rudimentary  sensitivity 
analysis  of  the  model's  solution  given  controlled  variation 
of  the  input  parameters.  This  chapter  explicitly  presents 
the  numerical  experiment  used  to  validate  the  integral 
model . 

The  general  approach  to  the  first  part  of  this 
numerical  experiment  is  as  follows:  A  modified  version  of 
the  SEARCH 4  program  is  used  to  create  hypothetical  masses 
on  three  airborne  filter  samples  given  fixed  input  param¬ 
eters.  Then,  these  sample  masses  are  used  as  input  to  the 
SEARCH4  program  for  analysis.  If  the  program  functions 
properly,  the  output  solution  vector,  consisting  of  three 
"A"  values  (fraction  of  mass  in  each  mass  group)  will  have 
a  minimum  standard  deviation  when  the  lognormal  distribution 
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parameters  being  checked  by  SEARCH4  match  the  lognormal 
parameters  used  to  create  the  samples. 

The  second  part  of  the  numerical  experiment  uses 
the  sample  data  created  in  the  first  part  as  a  base  case 
and  varies  input  parameters,  one  at  a  time,  to  study  the 
effect  on  the  output  solution  vector.  The  parameters  to  be 
varied  are  chosen  because  they  are  the  ones  that  are  most 
likely  to  be  different  in  an  actual  cloud  sampling  situa¬ 
tion  from  what  was  assumed  in  this  model.  These  variable 
parameters  include  the  falling  particle  density;  the  weapon 
yield;  the  total  mass  lofted  by  the  weapon,  and  the  follow¬ 
ing  sample  input  data:  sample  total  mass,  time  of  collec¬ 
tion,  and  collection  altitude. 

Further  details  concerning  the  numerical  experiment 
and  results  are  contained  in  the  following  sections  of  this 
chapter . 

Integral  Method  Numerical 
Validation  Experiment 

The  modifications  made  to  the  SEARCH4  program  in 
order  to  make  it  produce  hypothetical  filter  sample  data 
can  best  be  described  by  using  the  nine  step  program 
algorithm  presented  in  the  final  section  of  Chapter  II. 

The  procedure  involves  changing  steps  one,  two,  three,  and 
five  of  the  algorithm,  and  eliminating  steps  seven  through 
nine  completely.  Step  one  is  changed  so  that  the  need  to 
input  a  total  sample  mass  is  eliminated.  Step  two  is 
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changed  to  an  interactive  sequence  that  requests  a  median 
radius  and  log-slope  for  the  "known"  lognormal  particle 
number-size  distribution  to  be  used  throughout  the  numeri¬ 
cal  experiment.  Step  three  is  modified  so  that  the  dis¬ 
tribution  is  divided  into  50  equal  mass  groups.  Finally, 
step  five  is  altered  in  two  ways.  The  masses  contributed 
by  each  equal  mass  group  are  summed  for  each  different 
sample  time  and  altitude  and  stored  in  an  array.  Each  of 
these  values  represents  a  hypothetical  total  mass  on  an  air¬ 
borne  filter  sample  taken  at  a  given  altitude  and  a  given 
time.  Last  of  all,  these  newly  created  airborne  filter 
sample  masses  are  written  to  the  ANSWER  file  for  output. 

The  decision  to  use  50  equal  mass  groups  to  create 
the  hypothetical  filter  sample  data  is  based  on  numerical 
experimentation  fostered  by  the  following  reasoning.  Since 
in  reality,  the  stabilized  nuclear  cloud  contains  a  great 
number  of  particle  sizes,  50  equal  mass  groups  represented 
by  50  different  median  radii  more  validly  approximate  the 
actual  falling  cloud  than  only  a  few  equal  mass  groups 
with  only  a  few  median  radii.  Additionally,  one  of  the 
goals  of  this  validation  experiment  was  to  test  the  pro¬ 
gram's  ability  to  analyze  actual  sample  data.  It  was 
learned  through  numerical  experimentation  that  if  filter 
sample  masses  were  created  using  only  a  few  equal  mass 
groups,  thi„  procedure  not  only  failed  to  validly  approxi¬ 
mate  nature,  but  also  failed  to  challenge  the  program. 
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Consequently,  numerical  experimentation  was  used 
to  determine  the  optimum  number  of  equal  mass  groups  to 
use  in  creating  realistic  filter  sample  data.  Filter 
sample  masses  were  created  with  a  varying  number  of  equal 
mass  groups.  Sample  masses  created  with  3  to  15  equal  mass 
groups  varied  numerically  in  the  first  and  second  signifi¬ 
cant  figures.  Those  sample  masses  created  using  15  to  50 
equal  mass  groups  varied  in  the  second  and  third  signifi¬ 
cant  figure,  and  beyond  50  equal  mass  groups,  the  differ¬ 
ences  were  in  the  fourth  and  fifth  significant  figures. 
Therefore,  it  can  be  argued  that  50  equal  mass  groups  can 
be  used  to  adequately  represent  the  actual  falling  nuclear 
cloud  for  the  purposes  of  this  experiment. 

The  Defense  Land  Fallout  Information  Code  (DELFIC) 
default  spectrum  lognormal  parameters,  median  radius, 
rQ  =  0.204  micrometers,  and  log-slope,  B  =  ln(4)  (4:210) 

are  used  throughout  this  numerical  experiment.  This  choice 
is  made  for  no  reason  other  than  the  fact  that  the  DELFIC 
default  distribution  is  one  that  is  widely  used  in  fallout 
modeling.  The  other  input  is  arbitrarily  selected  as 
representative  of  a  typical  nuclear  cloud  sampling  project. 
Table  I  lists  the  input  data  for  the  modified  version  of 
SEARCH 4  and  Table  II  lists  the  output  filter  sample  data. 

The  calculated  filter  sample  masses  listed  in 
Table  II  are  considerably  larger  than  expected.  This  is 
because  they  represent  a  clean  sweep  of  all  the  mass  per 
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TABLE  I 

PROGRAM  INPUT  DATA  FOR  INTEGRAL  MODEL 
VALIDATION,  PART  ONE 


Yield 

(mt) 


Total  Mass  Lofted 
(kg) 


_ Sample _ 

Altitude  Time 

(m)  (min) 


i • 

13100 

182 

3.5 

1.06  x  10^  (*) 

15800 

152 

1 

16200 

184 

(*)  Note:  Mass  lofted  based  on  one-third  of  a  ton 
of  mass  lofted  per  ton  of  yield.  Additionally,  the  density 
of  this  material  is  assumed  to  be  2600  kilgrams  per  cubic 
meter . 


TABLE  II 

CALCULATED  HYPOTHETICAL  FILTER  SAMPLE  DATA 


Altitude  (m) 


13100 

15800 

16200 


Sample 


Time  (min) 


Total  Mass  (kg) 


41126 

47636 

42575 


£ 
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vertical  meter  in  the  nuclear  cloud  at  their  respective 
altitudes  and  at  their  respective  times  of  collection.  If 
the  assumptions  that  all  samples  were  taken  under  the  same 
conditions  and  the  aircraft  flew  through  the  geometric 
center  of  the  cloud  at  each  altitude  are  valid,  then  no 
correction  factors  are  needed  to  correct  for  different 
cloud  penetration  flight  paths  or  for  different  sampling 
air  flow  rates.  That  being  the  case,  the  magnitudes  of  the 
individual  sample  masses  are  not  important.  What  is  impor¬ 
tant  is  the  ratio  of  the  samples,  one  to  another. 

The  data  from  Tables  I  and  II  was  input  to  the 
SEARCH 4  program  and  a  quick  search  was  conducted  of  log¬ 
normal  distributions  with  log-slopes  of  1.5  to  5.0.  The 
first  search  revealed  no  solution  vectors  with  standard 
deviations  less  than  100.  In  this  case,  a  marginal  solu¬ 
tion  occurred  because  of  the  relatively  high  altitudes 
and  late  times  of  the  samples.  Most  of  the  large  particles 
had  already  fallen  from  the  cloud.  At  3  hours  after  stabil¬ 
ization,  for  a  3  megaton  burst,  particles  with  radii  of 
80  micrometers  are  hitting  the  ground  (6:67).  Therefore, 
the  majority  of  the  particles  larger  than  80  micrometers 
have  already  fallen  to  earth,  or  at  the  very  least,  are 
not  present  at  the  high  sample  altitudes.  Yet,  by  using 
only  three  equal  mass  groups  in  the  search,  the  model  is 
limited  to  working  with  group  median  radii  of  17,  65,  and 
250  micrometers.  It's  doubtful  that  any  250  micrometer 


particles  are  at  these  high  sample  altitudes  at  3  hours 
after  cloud  stabilization.  Consequently,  a  marginal  solu¬ 
tion  was  understandable. 

One  solution  to  the  large  particle  dilemma  is  to 
divide  the  mass  aloft  into  more  equal  mass  groups.  Then 
only  the  first  three  groups,  or  the  ones  with  the  smallest 
median  radii  are  used.  For  example,  in  this  experiment, 
the  mass  moment  of  the  lognormal  distribution  representing 
the  total  mass  aloft  was  divided  into  six  equal  mass 
groups.  Their  mean  radii  were,  9.6,  25.6,  48.6,  87,  166, 
and  443  micrometers.  Since  the  majority  of  the  particles 
with  the  three  largest  radii  have  fallen  well  below  the 
sampling  altitudes  during  the  3  hours  prior  to  sampling, 
the  50  percent  of  the  lofted  mass  contained  in  these 
groups  does  not  contribute  to  the  samples.  Therefore, 
since  only  the  first  three  groups  are  contributing  to  the 
samples,  and  they  contain  only  50  percent  of  the  mass 
aloft,  each  of  the  solution  vector  values  should  be  0.1667, 
or  one-sixth  of  the  mass  aloft. 

Before  any  program  output  is  presented,  an  explana¬ 
tion  of  the  method  of  presentation  is  necessary.  All  pro¬ 
gram  output  tabulated  throughout  this  study  is  presented 
in  the  following  format:  Each  potential  solution  distribu¬ 
tion  is  tabulated  in  terms  of  its  characteristic  parameters, 
median  radius  and  log-slope.  Additionally,  the  solution 
vector  of  "A"  values  and  a  standard  deviation  of  these  "A" 


values  from  their  algebraic  mean  are  presented  for  every 
potential  solution  distribution. 

Table  III  contains  the  output  from  a  search  con¬ 
ducted  with  the  program  SEARCH4  and  six  equal  mass  groups. 
Nine  potential  solutions  were  found  by  the  program.  How¬ 
ever,  the  solution  with  beta  =  4  and  the  median  radius 
=  0.2  is  the  only  one  that  has  two  solution  vector  elements 
closest  to  the  expected  value  of  0.1667. 

Two  additional  searches  were  conducted  in  an 
attempt  to  pinpoint  the  distribution  with  the  closest  fit. 
One  search  used  seven  equal  mass  groups  and  the  other  used 
eight.  Also,  both  of  these  searches  used  a  finer  mesh  of 
lognormal  median  radii,  thereby  checking  ten  times  as  many 
distributions  as  were  checked  by  the  previous  six  group 
search. 

Since  the  total  mass  aloft  is  now  being  divided 
into  seven  and  eight  equal  mass  groups,  the  expected  solu¬ 
tion  vector  values  are  0.143  and  0.125  respectively.  Also, 
by  dividing  the  mass  aloft  into  a  greater  number  of  groups, 
and  then  using  only  the  first  three  groups  for  the  sample 
analysis,  implicitly,  more  mass  is  assumed  to  have  fallen 
from  the  cloud  and  not  to  have  contributed  to  the  samples. 
For  example,  in  the  seven  group  analysis,  57  percent  of 
the  mass  is  assumed  to  be  on  the  ground,  or  at  least  below 
the  sampling  altitudes  prior  to  the  sampling  times.  Like¬ 
wise,  for  the  eight  group  analysis,  62  percent  of  the  mass 


TABLE  III 


SIX  EQUAL  MASS  GROUP  PARTICLE  DISTRIBUTION  SEARCH  RESULTS 


INPUT  SAMPLE  DATA 


SAMPLE 

#1: 

MASS  = 

41126 

ALTITUDE 

=  13100 

TIME  = 

182 

SAMPLE 

#2: 

MASS  = 

47636 

ALTITUDE 

=  15800 

TIME  = 

152 

SAMPLE 

#3: 

MASS  = 

42575 

ALTITUDE 

=  16200 

TIME  = 

184 

MEDIAN 

RADIUS 

(um) 

BETA 

STANDARD 

DEVIATION 

2.0 

2.9 

.  493  5933E-01 

SOLUTION  VECTOR 
. 2135982E+00 
. 1254121E+00 
. 1310807E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

1.0  3.1  . 1432091E-01 


SOLUTION  VECTOR 
•1442953E+00 
. 1344202E+00 
. 1160741E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

1.0  3.2  . 2219986E-01 

SOLUTION  VECTOR 
. 1822619E+00 
. 1421405E+00 
. 1457329E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

1.0  3.3  . 4534747E-01 


SOLUTION  VECTOR 
. 2292991E+00 
.1386594E+00 
. 186  718  7E+00 


TABLE  III--Continued 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.50  3.5  . 1805175E-01 


SOLUTION  VECTOR 
. 1453390E+00 
. 1644582E+00 
. 1283  762E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.30  3.8  . 436  3031E-01 


SOLUTION  VECTOR 
. I665152E+00 
. 158  712  7E+00 
. 2378811E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.20  4.0  . 5155210E-01 


SOLUTION  VECTOR 
. 1577227E+00 
. 1657446E+00 
. 2507538E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.11  4.3  . 6384798E-01 


SOLUTION  VECTOR 
. 1442388E+00 
. 1768401E+00 
. 2674626E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.090  4.4  . 6300539E-01 


SOLUTION  VECTOR 
. 1370042E+00 
.  1825519E+00 
.  2615282E+00 


is  assumed  to  be  below  the  sample  altitudes  at  sampling 
times.  These  assumptions  may  be  acceptable  at  later  times 
but  they  become  less  valid  as  the  time  between  cloud 
stabilization  and  sampling  is  decreased.  At  these  early 
times  following  stabilization  many  larger  particles  are 
still  aloft. 

In  the  seven  group  search,  the  number  of  potential 
solution  vectors  was  reduced  from  the  nine  listed  in 
Table  III  to  the  six  listed  in  Table  IV.  However,  the 
solution  vector  for  the  distribution  with  the  log-slope 
equal  to  four  had  a  slightly  lower  standard  deviation  than 
any  of  the  others  listed  in  the  table.  On  the  other  hand, 
the  individual  solution  vector  elements  for  this  distribu¬ 
tion  do  not  closely  approach  their  expected  values  of  0.143 
as  they  did  in  the  six  group  search.  Nonetheless,  based 
on  the  data  in  Tables  III  and  IV,  the  best  estimate  of  the 
actual  stabilized  cloud  distribution  at  this  point  is  that 
it  is  lognormal  with  a  log-slope  somewhere  between  3.9 
and  4.1,  and  a  median  radius  between  0.19  and  0.28  micro¬ 
meters  . 

The  last  search  conducted  is  the  eight  group  search. 
In  this  search,  the  number  of  potential  solution  vectors 
was  reduced  from  the  six  listed  in  Table  IV  to  the  four 
listed  in  Table  V.  This  time,  all  the  solution  vectors' 
elements  are  close  to  the  expected  value  of  0.125,  with  the 
exception  of  the  last  distribution.  Overall,  the  potential 


TABLE  IV 


RESULTS  FROM  THE  SEVEN  EQUAL  MASS  GROUP  SEARCH 


INPUT  SAMPLE  DATA 


SAMPLE 

#1: 

MASS  = 

41126 

ALTITUDE 

=  13100 

TIME  = 

182 

SAMPLE 

#2: 

MASS  = 

47636 

ALTITUDE 

=  15800 

TIME  = 

152 

SAMPLE 

#3: 

MASS  = 

42575 

ALTITUDE 

=  16200 

TIME  = 

184 

MEDIAN 

RADIUS 

(um) 

BETA 

STANDARD 

DEVIATION 

0.42 

3.6 

.8528  26  7E-02 

SOLUTION  VECTOR 
.1362968E+00 
. 136  786  5E+00 
. 1217763E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.34  3.7  . 7654989E-02 


SOLUTION  VECTOR 
. 1313907E+00 
. 1387913E+00 
.1234842E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.34  3.8  . 9797272E-02 


SOLUTION  VECTOR 
. 1598086E+00 
. 1563453E+00 
. 14 13  74  8E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.28  3.9  . 6429054E-02 


SOLUTION  VECTOR 
. 15  702  77E+0Q 
. 1594050E+00 
. 14  72728E+00 


TABLE  IV — Continued 


MEDIAN  RADIUS 
0.23 

(urn) 

BETA 

4.0 

STANDARD  DEVIATION 
. 6360630E-02 

SOLUTION  VECTOR 
. 1531412E+00 
. 16  29  241E+00 
. 1509903E+00 

MEDIAN  RADIUS 
0.19 

(um) 

BETA 

4.1 

STANDARD  DEVIATION 
. 747499E-02 

SOLUTION  VECTOR 
. 1507499E+00 
. 1656834E+00 
-1576448E+00 

TABLE  V 


# 


+ 


RESULTS  FROM  THE  EIGHT  EQUAL  MASS  GROUP  SEARCH 


INPUT  SAMPLE  DATA 


SAMPLE 

#1: 

MASS  =  41126 

ALTITUDE  =  13100 

TIME  =182 

SAMPLE 

#2: 

MASS  =  47636 

ALTITUDE  =  15800 

TIME  =152 

SAMPLE 

#3: 

MASS  =  42575 

ALTITUDE  =  16200 

TIME  =  184 

MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.27  3.9  . 7414105E-02 


SOLUTION  VECTOR 
. 1391024 E+ 00 
. 13  276  21E+00 
. 1243  23  7E+00 


MEDIAN  RADIUS  (tom)  BETA  STANDARD  DEVIATION 

0.22  4.0  . 4736288E-02 


SOLUTION  VECTOR 
. 1358764E+00 
. 1329348E+00 
. 1266077E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.18  4.1  . 3188375E-02 


SOLUTION  VECTOR 
. 1324259E+00 
. 1343090E+00 
. 1280913 E+ 00 


» 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.19  4.2  . 9649241E-02 


SOLUTION  VECTOR 
. 1696385E+00 
. 1531179E+00 
. 1700169E+00 


# 
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solutions  found  by  SEARCH4  during  this  eight  group  search 
echo  the  findings  of  the  seven  group  search.  Namely,  the 
particle  size  distribution  for  the  hypothetical  stabilized 
nuclear  cloud  that  the  calculated  filter  samples  came  from 
is  a  lognormal  distribution  with  a  log-slope  between  4.0 
and  4.1  and  a  median  radius  between  0.18  and  0.22  micro¬ 
meters  . 

In  summation,  the  integral  method  of  unfolding  the 
stabilized  nuclear  cloud  particle  size  distribution  from 
airborne  filter  sample  data  was  tested  in  this  section 
using  filter  sample  data  created  from  a  hypothetical 
nuclear  cloud.  This  hypothetical  nuclear  cloud's  particle 
size  distribution  was  the  DELFIC  default  particle  size 
distribution.  Since  the  SEARCH4  program  predicted  a  par¬ 
ticle  size  distribution  that  was  quite  close  to  the  "known" 
distribution,  used  to  create  the  filter  samples,  it  can  be 
argued  that  the  integral  model  is  a  valid  method  for  deter¬ 
mining  the  particle  size  distribution  of  a  stabilized  cloud, 
or  at  least  a  valid  method  for  determining  the  distribution 
log-slope . 

The  next  section  of  this  chapter  presents  a 
detailed  sensitivity  analysis  of  the  model  using  the  data 
in  Table  V  as  base  case  data. 


Integral  Method  Numerical 
Sensitivity  Analysis 

Many  of  the  measurements  made  during  the  atmos¬ 
pheric  nuclear  tests  produced  data  that  were  not  as  accu¬ 
rate  or  as  well  documented  as  that  which  is  needed  for  the 
integral  model.  For  example,  perhaps  the  actual  flight 
path  of  the  sampling  aircraft  did  not  pass  through  the 
geometric  center  of  the  nuclear  cloud.  In  most  cases,  no 
specific  notes  on  the  flight  path  were  made  aside  from  air¬ 
craft  altitude  and  time  in  the  cloud.  Also,  many  clouds 
from  megaton  sized  weapons  were  only  peripherally  sampled. 
Since  this  model  assumes  a  flight  path  through  the  geo¬ 
metric  center  of  the  nuclear  cloud  at  every  sampling  alti¬ 
tude,  a  difference  between  the  actual  flight  path  and  the 
assumed  flight  path  may  change  the  program  output  con¬ 
siderably  . 

Consequently,  the  purpose  of  this  part  of  the  inte¬ 
gral  model  testing  is  to  determine  how  sensitive  the  model 
is  to  small  inaccuracies  in  the  input  parameters.  The  cal¬ 
culated  filter  sample  masses  from  Table  II,  and  the  results 
of  the  eight  equal  mass  group  search  from  Table  V  are  used 
as  a  base  case  for  this  sensitivity  analysis.  Each  of  the 
program  input  parameters  is  varied,  one  at  a  time,  and  the 
effects  on  the  output  distribution  parameters  and  solution 
vectors  are  analyzed. 
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The  first  parameter  to  be  varied  is  the  particle 
density.  The  particle  density  used  in  calculating  the 
filter  sample  masses  found  in  Table  II  was  the  integral 
model  default  value,  2600  kilograms  per  cubic  meter.  In 
this  numerical  experiment,  the  density  input  to  SEARCH4  is 
varied  from  2300  kilograms  per  cubic  meter  to  2900  kilo¬ 
grams  per  cubic  meter.  The  resultant  distribution  param¬ 
eters  and  solution  vectors  for  the  search  conducted  using 
a  fallout  density  of  2300  kilograms  per  cubic  meter  are 
listed  in  Table  VI.  The  results  of  the  search  conducted 
using  a  fallout  density  of  2900  kilograms  per  cubic  meter 
are  listed  in  Table  VII.  When  the  assumed  particle  density 
is  varied  by  approximately  10  percent,  either  up  or  down, 
the  predicted  distribution  log-slope  and  median  radius 
remain  relatively  unchanged.  For  example,  from  Table  V, 
with  the  fallout  density  set  at  2600  kilograms  per  cubic 
meter,  the  optimum  distribution  is  one  with  a  log-slope 
between  3.9  and  4.2,  and  a  median  radius  between  0.18  and 
0.27  micrometers.  For  both  fallout  densities  of  2300  and 
2900  kilograms  per  cubic  meter,  the  optimum  distribution 
is  one  with  a  log-slope  between  3.6  and  4.2,  and  a  median 
radius  between  0.2  and  0.5  micrometers.  Therefore,  it  is 
evident  that  small  changes  in  the  fallout  particle  density 
do  not  significantly  effect  the  integral  irodel  output. 

This  is  as  expected  because  varying  the  particle  density 
by  10  percent  in  Equation  (3)  affects  all  falling  particles 
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TABLE  VI 


RESULTS  FROM  PARTICLE  DENSITY  VARIATION  NUMERICAL 
EXPERIMENT  FOR  DENSITY  =  2300  KILOGRAMS 
PER  CUBIC  METER 


BASE  CASE  SAMPLE  DATA 


SAMPLE 

#1: 

MASS  =  41126 

ALTITUDE  =  13100 

TIME  =182 

SAMPLE 

#2: 

MASS  =  47636 

ALTITUDE  =  15800 

TIME  =  152 

SAMPLE 

#3: 

MASS  =  42575 

ALTITUDE  =  16200 

TIME  =  184 

MEDIAN  RADIUS  (uni)  BETA  STANDARD  DEVIATION 

0.50  3.6  . 1295183E-01 


SOLUTION  VECTOR 
. 142908  5E+00 
. 12  24889E+00 
. 11889  54E+00 


MEDIAN  RADIUS  (um) 
0.50 


BETA  STANDARD  DEVIATION 

3.7  . 304348  2E-0I 


SOLUTION  VECTOR 
. 1663867E+00 
. 1491839E+00 
. 1072196E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.50  3.8  . 3852854E-01 


SOLUTION  VECTOR 
.21552  8  2E+00 
. 1387422E+00 
. 18  2728 1E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.30  3.9  . 3421003E-01 


SOLUTION  VECTOR 
. 1241103 E+ 00 
. 1663406E+00 
. 9860537E-01 


MEDIAN  RADIUS 

0.20 


(um) 


BETA 

4.0 


STANDARD  DEVIATION 
. 172  73  75E-01 


(um) 


(um) 


SOLUTION  VECTOR 
. 1784879E+00 
.  149  7153E+00 
.  1475411E+00 


BETA  STANDARD  DEVIATION 

4.1  . 4048359E-01 


SOLUTION  VECTOR 
.  2236078E+00 
.  1468642E+00 
. 162884  7E+00 


BETA  STANDARD  DEVIATION 

4.2  . 1022544E-01 


SOLUTION  VECTOR 
. 1655476 E+ 00 
. 15774 15E+00 
. 1452745E+00 


TABLE  VII 


RESULTS  FROM  PARTICLE  DENSITY  VARIATION  NUMERICAL 
EXPERIMENT  FOR  DENSITY  =  2900  KILOGRAMS 
PER  CUBIC  METER 


BASE 

CASE  SAMPLE  DATA 

SAMPLE 

#1: 

MASS  = 

41126 

ALTITUDE  =  13100 

TIME  = 

182 

SAMPLE 

#2: 

MASS  = 

47636 

ALTITUDE  =  15800 

TIME  = 

152 

SAMPLE 

#3: 

MASS  = 

42575 

ALTITUDE  =  16200 

TIME  = 

184 

MEDIAN 

RADIUS 

(um) 

BETA  STANDARD 

DEVIATION 

0.50 

3.6  . 1956538E-01 

SOLUTION  VECTOR 

• 

1578669E+00 

• 

1342443 E+ 00 

• 

1190391E+00 

MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.50  3.7  . 3570383E-01 


SOLUTION  VECTOR 
. 1994141E+00 
. 1337939E+00 
.  1909921E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.30  3.9  . 1256950E-01 


SOLUTION  VECTOR 
.  1619  730E+00 
.  1467839E+00 
.  1370307E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.20  4.1  . 9405460E-02 


SOLUTION  VECTOR 
. 1501645E+00 
. 1543116E+00 
. 1363482E+00 


equally,  and  therefore  would  have  little  effect  on  the  out¬ 
put  from  SEARCH4 . 

The  second  parameter  to  be  varied  is  the  yield. 
Since  the  cloud  rise  portion  of  the  model  is  heavily  depen¬ 
dent  upon  the  weapon  yield  in  order  to  position  the  par¬ 
ticle  groups  at  their  initial  stabilized  altitudes,  a 
difference  between  the  actual  yield  and  the  assumed  yield 
could  affect  the  program's  output. 

For  this  numerical  experiment,  the  yield  used  to 
create  the  filter  samples  found  in  Table  II  was  varied  from 
3  megatons  to  4  megatons.  The  resultant  distribution 
parameters  and  solution  vectors  are  listed  in  Tables  VIII 
and  IX.  Again,  no  significant  changes  in  the  SEARCH4  out¬ 
put  are  discernible  as  the  yield  is  varied  up  to  15  per¬ 
cent  in  either  direction.  The  SEARCH4  program  still 
selects  an  optimum  distribution  with  a  log-slope  between 
3.7  and  4.4,  and  a  median  radius  between  0.1  and  0.5  micro¬ 
meters  . 

Regardless  of  a  slight  yield  variation,  stable 
program  output  is  understandable  for  the  following  reasons. 
In  this  case,  an  increase  or  decrease  in  yield  of  half  a 
magaton  only  changes  the  particles'  initial  altitudes  by 
about  500  meters  for  particles  in  the  10  to  150  micrometer 
size  range.  This  small  initial  altitude  change  produces 
only  a  minor  difference  in  the  mass  per  vertical  meter  at 
the  sample  altitudes.  Additionally,  given  one  set  of 


TABLE  VIII 


RESULTS  FROM  WEAPON  YIELD  VARIATION  NUMERICAL  EXPERIMENT 

FOR  YIELD  =  3  MEGATONS 


BASE  CASE  SAMPLE  DATA 


SAMPLE 

#1: 

MASS  = 

41126 

ALTITUDE 

=  13100 

TIME  = 

182 

SAMPLE 

#2: 

MASS  = 

47636 

ALTITUDE 

=  15800 

TIME  = 

152 

SAMPLE 

#3: 

MASS  = 

42575 

ALTITUDE 

=  16200 

TIME  = 

184 

MEDIAN 

RADIUS 

(um) 

BETA 

STANDARD 

DEVIATION 

0.2  4.0  . 2562264E-01 


SOLUTION  VECTOR 
. 1302901E+0Q 
. 1353253E+00 
. 8864276E-01 


MEDIAN  RADIUS 
0.2 


(tun)  BETA  STANDARD  DEVIATION 

4.2  . 4721363E-01 


SOLUTION  VECTOR 
. 2074971E+00 
.  1237366E+00 
. 1278609E+00 


I 


MEDIAN  RADIUS 

0.11 


(um)  BETA  STANDARD  DEVIATION 

4.3  . 224  7532E-01 


SOLUTION  VECTOR 
. 1204685E+00 
. 13  919  76E+00 
. 9444473E-01 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.09  4.4  . 1906333E-01 


SOLUTION  VECTOR 
. 1196479E+00 
. 1364339E+00 
. 9839454E-01 


TABLE  IX 


RESULTS  FROM  WEAPON  YIELD  VARIATION  NUMERICAL  EXPERIMENT 

FOR  YIELD  =  4  MEGATONS 


BASE  CASE  SAMPLE  DATA 


SAMPLE 

#1: 

MASS  = 

41126 

ALTITUDE 

=  13100 

TIME  = 

182 

SAMPLE 

#2: 

MASS  = 

47636 

ALTITUDE 

=  15800 

TIME  = 

152 

SAMPLE 

#3: 

MASS  = 

42575 

ALTITUDE 

=  16200 

TIME  = 

184 

MEDIAN 

RADIUS 

(um) 

BETA 

STANDARD 

DEVIATION 

0.5 

3.7 

.  10906  53E-01 

SOLUTION  VECTOR 
. 1683339E+00 
. 1532609E+00 
. 1744525E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.3  3.9  . 9227757E-02 


SOLUTION  VECTOR 
. 1459867E+00 
. 1428191E+00 
. 1601487E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.3  4.0  . 6607308E-01 


SOLUTION  VECTOR 
. 18  2468  7E+00 
. 1536467E+00 
.2797444E+00 


TABLE  IX--Continued 


MEDIAN  RADIUS  (um) 
0.2 


BETA  STANDARD  DEVIATION 

4.1  . 2396100E-01 


SOLUTION  VECTOR 
.  1536515E+00 
.  1277581E+00 
•1756267E+00 


MEDIAN  RADIUS  (um) 
0.2 


BETA  STANDARD  DEVIATION 

4.2  .5306  501E-01 


SOLUTION  VECTOR 
. 1681936E+00 
.1626484E+00 
.2572067E+00 


particle  size  groups  started  at  an  altitude  of  16000  meters 
and  a  like  set  of  groups  started  at  an  altitude  of  16000 
meters  plus  or  minus  500  meters,  the  mass  per  vertical 
meter  for  each  set  of  groups  is  in  the  same  ratio  at  any 
given  sample  altitude,  at  any  given  time. 

The  third  input  parameter  to  be  studied  is  the 
total  mass  of  material  lofted  by  the  weapon.  For  each  of 
the  numerical  experiments  in  this  study,  the  empirical 
value  of  one-third  of  a  ton  of  material  lofted  per  ton  of 
yield  is  used.  However,  the  total  amount  of  mass  lofted 
makes  no  difference  in  the  outcome  of  these  calculations. 
This  is  intuitively  obvious  from  Equations  (12),  (13), 

(14) ,  and  (15)  .  Mathematically,  the  total  mass  lofted, 

Mfc,  is  only  a  scalar  constant  that  is  multiplied  times  the 
"F"  matrix  in  Equation  (15)  and  has  no  effect  on  the  deter¬ 
mination  of  the  solution  vector  elements  other  than  to 
increase  or  decrease  them  by  the  same  amount. 

The  fourth  input  parameter  to  be  varied  is  the 
mass  on  the  filter  samples.  As  stated  earlier,  the  inte¬ 
gral  model  assumes  that  all  samples  were  taken  under  a  set 
of  standard  sampling  conditions.  These  conditions  include, 
a  flight  path  through  the  geometric  center  of  the  cloud 
at  every  sample  altitude,  and  an  equal  amount  of  airflow 
through  the  filter  sampling  device  at  all  sampling  alti¬ 
tudes.  Since  the  nuclear  cloud  is  assumed  to  be  normally 
distributed  in  the  horizontal  direction  as  well  as  in  the 


* 


vertical  direction,  any  sample  technique  other  than  a 
single  flight  through  the  geometric  center  of  the  cloud 
will  require  a  correction  factor  to  correct  the  sample  mass 
to  the  standard  sampling  conditions. 

In  order  to  study  the  effect  of  nonstandard  sam¬ 
pling  conditions  such  as  cloud  edge  sampling  or  multiple 
flights  through  the  cloud  at  the  same  altitude  with  the 
same  filter,  the  base  case  sample  mass  data  is  altered  in 
the  following  two  ways.  First  the  effect  of  one  nonstan¬ 
dard  filter  sample  is  studied  by  decreasing  one  of  the 
sample  masses  by  3  percent  of  its  true  value  and  running 
this  data  through  SEARCH4 .  Then,  the  same  sample  is 
increased  by  10  percent  of  its  true  value  and  again  run 
through  SEARCH 4 .  The  resultant  distribution  parameters 
and  solution  vectors  are  listed  in  Tables  X  and  XI. 

The  results  in  Table  X  indicate  that  a  3  percent 
error  in  one  of  the  sample  masses  makes  very  little  differ¬ 
ence  in  the  optimum  particle  size  distribution  that  the 
SEARCH 4  program  selects.  Moreover,  the  data  in  Table  XI 
clearly  indicates  that  as  much  as  10  percent  error  in  one 
of  the  filter  sample  masses  has  only  minor  impact  on  the 
optimum  particle  size  distribution  selected  by  SEARCH4 . 

The  program  continues  to  select  an  optimum  particle  size 
distribution  that  is  lognormal  with  a  log-slope  between 
3.7  and  4.0,  and  a  median  radius  between  0.3  and  0.5  micro¬ 
meters. 
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TABLE  X 

RESULTS  FROM  VARIATION  OF  ONE  SAMPLE  MASS 
DOWN  3  PERCENT  FROM  41126  TO  39892 


BASE  CASE  SAMPLE  DATA 


SAMPLE 

SAMPLE 

SAMPLE 

#1: 

#2: 

#3: 

MASS  = 
MASS  = 
MASS  = 

41126 

47636 

42575 

ALTITUDE 

ALTITUDE 

ALTITUDE 

=  13100 
=  15800 
=  16200 

TIME  = 
TIME  = 
TIME  = 

182 

152 

184 

MEDIAN 

RADIUS 

(um) 

BETA 

STANDARD 

DEVIATION 

0.3  3.7  . 8829385E-01 


SOLUTION  VECTOR 
. 1539860E+00 
. 1656102E-01 
•1813110E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.2  3.9  . 784  7556E-01 

SOLUTION  VECTOR 
. 1388279E+00 
. 2916  525E-01 
. 1812376E+00 


TABLE  XI 


RESULTS  FROM  VARIATION  OF  ONE  SAMPLE  MASS 
UP  10  PERCENT  FROM  41126  TO  45239 


BASE  CASE  SAMPLE  DATA 


SAMPLE 

#1: 

MASS  = 

41126 

ALTITUDE 

=  13100 

TIME  = 

182 

SAMPLE 

#2: 

MASS  = 

47636 

ALTITUDE 

=  15800 

TIME  = 

152 

SAMPLE 

#3: 

MASS  = 

42575 

ALTITUDE 

=  16200 

TIME  = 

184 

MEDIAN 

RADIUS 

(um) 

BETA 

STANDARD 

DEVIATION 

0.5  3.7  . 3525377E+00 

SOLUTION  VECTOR 
.  2548912E+00 
.  9350324E-02 
.  7045112E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

•  0.3  4.0  . 6636146E+00 

SOLUTION  VECTOR 
.  2424191E+00 
. 4390633 E- 01 
. 12  796  47E+01 


For  the  next  part  of  the  mass  sensitivity  analysis, 
two  of  the  sample  masses  are  increased  by  10  percent  of 
their  true  values.  Table  XII  contains  the  results  of  this 
analysis.  It  is  clear  that  if  two  sample  masses  are  inaccu 
rate  by  10  percent  in  the  same  direction,  the  resultant 
solution  vector  is  radically  changed  to  include  one  ele¬ 
ment  that  is  five  or  six  order  of  magnitude  above  the  other 
two  solution  vector  elements.  However,  it  is  most  impor¬ 
tant  to  note  that  the  optimum  distribution  selected  by 
SEARCH 4  remains  relatively  unchanged  with  a  log-slope 
between  3.7  and  4.0,  and  a  median  radius  between  0.5  and 
0.9  micrometers. 

The  final  mass  experiment  involves  varying  one  of 
the  sample  masses  by  plus  10  percent  of  its  true  value, 
and  varying  a  second  sample  mass  by  minus  10  percent  of 
its  true  value.  When  this  is  done,  no  good  solution  is 
found  by  SEARCH4.  Therefore,  it  is  clearly  evident  that 
relatively  small  inaccuracies  in  the  input  sample  masses 
can  greatly  influence  the  integral  model  output.  This 
extreme  sensitivity  to  individual  sample  mass  differences 
is  understandable  since  the  whole  basis  for  the  optimum 
distribution  selection  in  the  integral  model  is  that  the 
total  mass  on  the  filter  is  equal  to  the  sum  of  the  masses 
contributed  by  the  equal  mass  groups  being  used  to  repre¬ 
sent  the  falling  nuclear  cloud.  Therefore,  if  the  total 
mass  on  the  filter  is  not  accurately  acquired,  the  sum  of 


. 
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TABLE  XII 


i* 

V  RESULTS  FROM  VARIATION  OF  TWO  SAMPLE  MASSES — 
£  SAMPLE  #1:  UP  10  PERCENT  FROM  41126  to  45239; 
U  SAMPLE  #2:  UP  10  PERCENT  FROM  47636  to  52400 


the  mass  group  contributions,  calculated  by  SEARCH4 ,  will 
be  considerably  different  from  the  total  sample  mass  and  no 
optimum  distribution  will  be  found. 

The  fifth  input  parameter  to  be  varied  is  the 
sample  time  of  collection.  In  the  first  experiment,  one 
sample  time  is  varied  by  minus  10  percent  of  its  true 
value.  Next,  two  sample  times  are  varied  by  minus  10  per¬ 
cent  of  their  true  values.  The  resultant  SEARCH4  output 
is  contained  in  Tables  XIII  and  XIV.  It  is  apparent  that 
small  inaccuracies  in  the  sample  times  of  collection  make 
little  impact  on  the  SEAJRCH4  program  output.  This  is 
especially  true  at  late  times  when  most  of  the  larger  par¬ 
ticles  have  fallen  from  the  cloud  and  only  the  slowly  fall¬ 
ing  smaller  particles  are  still  airborne. 

The  final  input  parameter  to  be  varied  is  the 
sample  altitude.  First,  one  sample  altitude  is  increased 
by  5  percent  of  its  original  value.  Second,  two  sample 
altitudes  are  increased  by  5  percent  of  their  original 
values.  The  results  of  this  experiment  are  contained  in 
Tables  XV  and  XVI. 

The  Table  XV  data  indicates  that  one  sample  alti¬ 
tude  error  of  5  percent  or  less  does  not  significantly 
affect  the  program  output.  However,  from  the  data  in 
Table  XVI,  it  is  evident  that  if  two  sample  altitudes  are 
in  error  by  as  little  as  5  percent,  the  program  output  is 
distinctly  altered.  One  of  the  solution  vector  elements 
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TABLE  XIII 

RESULTS  FROM  VARIATION  OF  ONE  SAMPLE  TIME 
DOWN  10  PERCENT  FROM  182  TO  164 


BASE  CASE  SAMPLE  DATA 


SAMPLE 

#1: 

MASS  =  41126 

ALTITUDE  =  13100 

TIME  =182 

SAMPLE 

#2: 

MASS  =  47636 

ALTITUDE  =  15800 

TIME  =152 

SAMPLE 

#3: 

MASS  =  42575 

ALTITUDE  =  16200 

TIME  =  184 

MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.5  3.6  . 3232540E-01 


% 


SOLUTION  VECTOR 
. 1992990E+00 
. 1520953E+00 
. 13  74395E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.3  3.8  . 30 76463 E- 01 


SOLUTION  VECTOR 
. 15648  78E+00 
. 1750403E+00 
.  1149581E+00 


MEDIAN  RADIUS  (um) 
0.2 


BETA  STANDARD  DEVIATION 

4.0  . 3042188E-01 


SOLUTION  VECTOR 
. 1466314E+00 
. 18  290  93E+00 
. 1224689E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.11  4.3  . 3597496E-01 


SOLUTION  VECTOR 
.  131884  5E+00 
.194 7588 E+ 00 
.1330278E+00 
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TABLE  XIV 


RESULTS  FROM  VARIATION  OF  TWO  SAMPLE  TIMES — 
SAMPLE  #1:  DOWN  10  PERCENT  FROM  182  TO  164; 
SAMPLE  #2:  DOWN  10  PERCENT  FROM  152  TO  137 


BASE  CASE  SAMPLE  DATA 


SAMPLE 

#1: 

MASS  = 

41126 

ALTITUDE 

=  13100 

TIME  = 

182 

SAMPLE 

#2: 

MASS  - 

47636 

ALTITUDE 

=  15800 

TIME  = 

152 

SAMPLE 

#3: 

MASS  = 

42575 

ALTITUDE 

=  16200 

TIME  = 

184 

MEDIAN 

RADIUS 

(um) 

BETA 

STANDARD 

DEVIATION 

0.3  3.8  . 2269939E+00 


SOLUTION  VECTOR 
. 234433  5E+00 
. 5204663E-01 
. 5032817E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.2  4.0  . 233543  2E+00 


SOLUTION  VECTOR 
.2313492E+00 
. 5556365E-01 
. 518  2256E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.11  4.3  .2443441E+00 


SOLUTION  VECTOR 
. 2260754E+00 
. 6164088E-01 
.54  2396  7E+00 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.09  4.4  . 2411485E+00 


SOLUTION  VECTOR 
. 2245084E+00 
. 6268384E-01 
. 5370646E+00 


TABLE  XV 


RESULTS  FROM  VARIATION  OF  ONE  SAMPLE  ALTITUDE 


UP 

5 

PERCENT 

FROM  13100 

TO  13755 

BASE  CASE  SAMPLE  ! 

DATA 

SAMPLE 

#1: 

MASS 

41126 

ALTITUDE 

=  13100 

TIME  = 

182 

SAMPLE 

#2: 

MASS 

= 

47636 

ALTITUDE 

=  15800 

TIME  = 

152 

SAMPLE 

#3: 

MASS 

= 

42575 

ALTITUDE 

=  16200 

TIME  = 

184 

MEDIAN 

RADIUS 

(um) 

BETA 

STANDARD 

DEVIATION 

0.3  3.8  . 7183170E+0  0 


SOLUTION  VECTOR 
. 2  5  72  775E+00 
. 1056209E-01 
. 1359598E+01 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.2  4.0  . 6974788E+00 


SOLUTION  VECTOR 
. 2501462E+00 
. 2239275E-01 
. 1328128 E+ 01 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.11  4.3  . 6921787E+00 


SOLUTION  VECTOR 
. 2406420E+00 
. 36  73403E-01 
.13  24  500E+01 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.09  4.4  . 6669652E+00 


SOLUTION  VECTOR 
. 2381173E+00 
.39  738  70E-01 
. 1281299E+01 


TABLE  XVI 


RESULTS  FROM  VARIATION  OF  TWO  SAMPLE  ALTITUDES — 
SAMPLE  #1:  UP  5  PERCENT  FROM  13100  TO  13755; 
SAMPLE  #2:  UP  5  PERCENT  FROM  15800  TO  16590 


BASE 

CASE  SAMPLE  DATA 

SAMPLE 

#1: 

MASS  = 

41126 

ALTITUDE  = 

13100 

TIME  = 

182 

SAMPLE 

#2: 

MASS  = 

47636 

ALTITUDE  = 

15800 

TIME  = 

152 

SAMPLE 

#3: 

MASS  = 

42575 

ALTITUDE  = 

16200 

TIME  = 

184 

MEDIAN 

RADIUS 

(urn) 

BETA  STANDARD 

DEVIATION 

in 

• 

o 

3.8 

. 19126 12E+06 

SOLUTION  VECTOR 

.  2702619E+00 
.2028975E+00 
.  3312743E+06 

is  increased  by  six  orders  of  magnitude.  Nonetheless,  it 
is  important  to  note  that  the  SEARCH4  program  continues  to 
select  an  optimum  particle  size  distribution  with  log- 
slope  between  3.8  and  4.4,  and  median  radius  between  0.1 
and  0.5  micrometers. 

In  summary,  a  sensitivity  analysis  of  the  integral 
model  program,  SEARCH4,  was  conducted  in  the  following 
manner.  The  input  parameters  that  were  most  likely  to  be 
in  error  in  an  actual  nuclear  cloud  sampling  scenario  were 
varied,  one  at  a  time,  in  order  that  the  effect  on  the  pro¬ 
gram  output  be  evaluated.  Overall,  this  analysis  indi¬ 
cated  that  the  integral  model  program  is  extremely  sensi¬ 
tive  to  errors  in  the  airborne  sample  mass  collection  pro¬ 
cedure  and  to  sample  altitude  errors.  If  corrections  are 


not  made  for  nonstandard  sampling  procedures,  or  if  some  of 
the  sample  altitudes  are  not  accurate,  the  program  does  not 
definitively  pinpoint  the  median  radius  of  the  stabilized 
cloud  particle  size  distribution.  However,  it  is  impor¬ 
tant  to  note  that  in  all  cases,  SEARCH4  does  not  fail  to 
indicate  the  true  log-slope  of  the  stabilized  nuclear  cloud 
particle  distribution. 

This  chapter  presented  the  results  of  the  integral 
model  testing  and  evaluation,  given  hypothetical  airborne 
filter  sample  data  of  known  origin.  It  demonstrated  that 
the  integral  model  of  airborne  filter  sample  data  analysis 
is  a  viable  method  for  unfolding  the  stabilized  nuclear 
cloud  particle  size  distribution  from  the  total  mass  col¬ 
lected  on  airborne  filters.  The  next  chapter  presents  a 
study  of  actual  airborne  filter  sample  data  from  a  nuclear 
shot  in  the  Pacific  test  range.  Much  of  the  analysis  con¬ 
tained  in  this  chapter  will  apply  to  the  study  conducted 
in  the  following  chapter. 


IV.  Shot  ZUNI  Filter  Sample  Data  Analysis 
Background 

Actual  airborne  sample  data  from  United  States 
atmospheric  nuclear  test  shots  conducted  during  the  1950s 
and  early  1960s  is  not  easy  to  find.  One  reason  for  the 
scarcity  of  data  is  that  only  a  small  number  of  test  shot 
nuclear  clouds  were  sampled  by  aircraft.  In  most  cases, 
an  extensive  study  was  made  of  the  fallout  particulate 
material  and  radiation  exposure  rates  on  the  ground  down 
wind  from  the  nuclear  detonation,  but  cloud  samples  at  dif¬ 
ferent  altitudes  were  seldom  taken.  Another  reason  for  the 
difficulty  in  obtaining  airborne  sample  data  is  that  little 
research  requiring  nuclear  cloud  sample  data  has  been  done 
since  Nathans  published  his  work  in  1970  (19:360-371). 
Consequently,  most  of  what  little  data  exists  is  probably 
contained  in  the  classified  archives  at  the  national  labora 
tories  and  has  been  virtually  forgotten.  In  most  cases, 
the  type  of  data  that  is  required  by  this  study  is  not 
classified.  Such  is  the  case  for  the  shot  ZUNI  nuclear 
cloud  sample  data  cited  in  this  chapter. 

Nathans  listed  four  nuclear  test  shots  that  were 


sampled  by  aircraft  and  rockets.  These  shots  were  Castle- 
BRAVO  (1954),  Redwing-ZUNI  (1956),  Castle-KOON  (1954),  and 
Operation  Sunbeam-JOHNNY  BOY  (1962)  (19:362).  However, 


the  search  for  information  on  these  four  shots  yielded  only 
a  small  amount  of  airborne  sample  data  from  ZUNI.  In  the 
first  part  of  this  chapter,  the  ZUNI  airborne  sample  data 
is  presented  along  with  the  results  of  the  SEARCH4  computer 
analysis  of  this  data.  The  second  part  of  the  chapter  is 
devoted  to  an  analysis  of  the  results  of  the  SEARCH4  out¬ 
put  and  a  hypothesized  particle  size  distribution  for  the 
ZUNI  stabilized  cloud  based  on  those  results. 

ZUNI  Airborne  Sample  Data 
and  SEARCH 4  Results 

The  ZUNI  airborne  sample  data  is  contained  in 
Table  XVII  (Ref  5) .  There  are  some  unknown  facts  concern¬ 
ing  the  actual  method  used  for  this  sampling  that  lead  to 
uncertainties  in  the  input  parameters.  First  of  all,  this 
data  was  not  extracted  from  airborne  filter  samples. 

Rather,  it  is  believed  to  have  been  gathered  by  three  gas 
sampling  devices,  each  mounted  on  a  different  B-57B  air¬ 
craft.  Second,  aside  from  the  sample  altitude,  the  actual 
flight  paths  of  the  aircraft  with  respect  to  the  cloud  are 
unknown.  For  example,  it  is  not  known  if  the  aircraft  flew 
through  the  geometric  center  of  the  cloud  at  all  altitudes 
or  if  the  aircraft  sampled  only  the  cloud  edges.  Addi¬ 
tionally,  the  times  of  collection  suggest  that  different 
amounts  of  gas  were  collected  at  each  altitude.  Finally, 
the  accuracy  of  the  total  mass  per  sample  is  suspect  since 
it  is  rounded  to  the  nearest  whole  number.  All  of  these 


TABLE  XVII 


REDWING 

SER-lES,  SHOT 

ZUNI  AIRBORNE  SAMPLE 

DATA 

Flight 

Sample 

Altitude 

(ft/m) 

Sample 

Collection  Time 
(min) 

Sample 
Total  Mass 
(oz) 

TIGER  RED  2 

43000/13106 

171-193  (182)* 

9 

HOTSHOT  2 

51700/15758 

138-165  (151.5)* 

11 

KASSADY  1 

53300/16246 

169-199  (184)* 

5 

*  Note:  Mean  sample  time  of  collection. 


uncertainties  greatly  enhance  the  potential  for  the  exis¬ 
tence  of  nonstandard  sampling  conditions  like  those  pre¬ 
sented  in  Chapter  III.  But,  since  this  is  the  only  actual 
nuclear  cloud  sample  data  that  was  uncovered  during  the 
data  search,  it  constitutes  the  best  available  data  and 
merits  a  complete  analysis. 

In  order  to  properly  analyze  the  ZUNI  data  and 
better  interpret  the  analytical  results,  some  additional 
background  information  for  the  shot  is  necessary.  The 
device  was  detonated  on  a  platform  at  a  height  of  9  feet 
over  coral  soil  at  the  Bikini  atoll  on  28  May  1956  at  0556 
hours.  It  produced  a  yield  of  3.5  megatons.  The  nuclear 
cloud  stabilized  in  approximately  5  minutes.  Its  top 
reached  an  altitude  of  79000  feet  (24100  meters)  and  its 
bottom  stabilized  at  approximately  49000  feet  (14900 
meters) .  At  one  hour  following  the  detonation,  the 


tropopause  was  at  an  altitude  of  51000  feet  (15500  meters) 
(10:111,114).  Figure  3  graphically  portrays  the  ZUNI  cloud 
sampling  scenario  and  brings  into  perspective  the  sample 
locations  relative  to  the  nuclear  cloud  position  at 
sampling  time. 

The  ZUNI  airborne  sample  data  contained  in  Table 
XVIII  was  input  to  the  SEARCH4  program.  For  the  first  run, 
only  three  equal  mass  groups  were  used  for  the  calcula¬ 
tions.  The  results  from  this  run  are  contained  in  Table 
XIX.  Additional  SEARCH4  runs  were  made  with  four  through 
eight  equal  mass  groups.  The  results  from  these  additional 
runs  are  contained  in  Tables  XX  through  XXIV. 

One  change  is  made  in  order  to  accommodate  the  dif¬ 
ference  in  total  mass  per  vertical  meter  at  the  sampling 
altitudes  versus  the  very  small  fraction  of  that  total 
mass  collected  by  the  sampling  devices.  This  change  is 
necessary  to  produce  solution  vectors  that  are  directly 
comparable  to  the  ones  in  the  Chapter  III  sensitivity 
analysis.  In  that  sensitivity  analysis,  the  input  sample 
masses  represented  the  total  mass  per  vertical  meter  at 
the  sampling  altitudes.  The  ZUNI  samples,  listed  in  Table 
XVIII,  represent  only  a  very  small  fraction  of  the  total 
mass  per  vertical  meter  at  the  sampling  altitudes.  In 
fact,  they  are  five  orders  of  magnitude  less  than  the 
samples  created  for  the  Chapter  III  model  validation  study 
(listed  in  Table  II).  As  determined  during  the  model 
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TABLE  XVIII 


INPUT  DATA  FILE  "FILTER"  FOR  ZUNI  SAMPLE  ANALYSIS 


YIELD  =  3500  (kt)  NUMBER  OF  SAMPLES  =3 

DISTRIBUTION  MOMENT  =  3  MASS  ALOFT  =  1.06  x  104  (kg) 


(kq) 

(m) 

(min) 

SAMPLE 

#1: 

MASS  =  .255 

ALTITUDE  =  13106 

TIME  =  182.0 

SAMPLE 

#2: 

MASS  =  .3118 

ALTUTUDE  =  15758 

TIME  =  151.5 

SAMPLE 

#3: 

MASS  =  .142 

ALTUTUDE  =  16246 

TIME  =  184.0 

TABLE  XIX 

RESULTS  OF  ZUNI  DATA  SEARCH  USING  THREE  EQUAL  MASS  GROUPS 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

1.0  2.9  . 9479308E+06 


SOLUTION  VECTOR 
•3316389E-01 
. 1456210E+00 
.1641864E+07 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.05  4.8  . 1725676E+06 


SOLUTION  VECTOR 
.9411089E-01 
. 2988994E+06 
.  7072789E+01 


TABLE  XX 


RESULTS  OF  ZUNI  DATA  SEARCH  USING 
FOUR  EQUAL  MASS  GROUPS 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

5.0  2.5  . 6225378E+06 


SOLUTION  VECTOR 
. 1043  578E+00 
. 423  56  24E+00 
. 10  78  26  7E+0  7 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

4.0  2.6  . 9036186E+06 


SOLUTION  VECTOR 
. 9877644E-01 
. 4135901E+00 
.1565114E+07 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

3.0  2.7  . 12811 78E+06 


SOLUTION  VECTOR 
. 8797728E-01 
.2724262E+00 
. 22 1906  7E+06 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

1.0  3.2  . 8869613E+06 


SOLUTION  VECTOR 
. 7862210E-01 
.233569  5E+00 
. 1536262 E+ 07 


TABLE  XX--Continued 


MEDIAN  RADIUS  (um) 
0.5 


BETA  STANDARD  DEVIATION 

3.5  . 4926050E+06 


SOLUTION  VECTOR 
. 7075585E-01 
.  1847197E+00 
. 8532169E+06 


MEDIAN  RADIUS  (um) 
0.2 


BETA  STANDARD  DEVIATION 

3.9  . 1439821E+06 


SOLUTION  VECTOR 
. 5647343E-01 
.  1530068E+00 
.  2493845E+06 


TABLE  XXI 

RESULTS  OF  ZUNI  DATA  SEARCH  USING 
#  FIVE  EQUAL  MASS  GROUPS 


MEDIAN  RADIUS  (um) 

BETA 

STANDARD  DEVIATION 

00 

• 

o 

2.4 

.  1323250E+06 

SOLUTION  VECTOR 
.1451400E+00 
.2569704E+00 
.2291939E+06 

MEDIAN 

0 

RADIUS 

.06 

(um) 

BETA 

4.7 

STANDARD  DEVIATION 
. 1111025E+06 

• 

SOLUTION  VECTOR 
. 6348051E-01 
. 1655484E+00 
. 1924353E+06 

MEDIAN 

0 

RADIUS 

.05 

(um) 

BETA 

4.8 

STANDARD  DEVIATION 
. 1 7536  79E+06 

• 

SOLUTION  VECTOR 
. 6393483E-01 
. 1671012E+00 
. 303  746  2E+06 

MEDIAN 

0 

RADIUS 

.035 

(um) 

BETA 

5.0 

STANDARD  DEVIATION 
. 5701899E+06 

SOLUTION  VECTOR 
. 65366  77E-01 
.  1719004E-1-00 
. 9875980E+06 


TABLE  XXII 


RESULTS  OF  ZUNI  DATA  SEARCH  USING 
SIX  EQUAL  MASS  GROUPS 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

1.5  3.3  . 938946  7E+06 


SOLUTION  VECTOR 
. 1006949E+00 
. 7017220E+00 
. 1626304 E+ 07 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.5  3.8  . 121396  7E+06 


SOLUTION  VECTOR 
. 8438121E-01 
. 3214669E+00 
. 2102655E+06 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.07  4.8  . 2197830E+06 


SOLUTION  VECTOR 
. 7509261E-01 
. 2227036E+00 
. 3806755E+06 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.06  4.9  . 8750562E+06 


SOLUTION  VECTOR 
. 773  752  5E-01 
. 2541977E+00 
. 151564  2E+0  7 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

0.047  5.0  . 1244411E+06 


SOLUTION  VECTOR 
.715514  7E-0 1 
. 1963516E+00 
. 2155385E+06 


♦ 


TABLE  XXIII 

RESULTS  OF  ZUNI  DATA  SEARCH  USING 
SEVEN  EQUAL  MASS  GROUPS 


^  MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

5.0  2.8  . 4525707E+06 

SOLUTION  VECTOR 
. 1549908E+00 
. 1833116E+00 

«  . 7838756E+06 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

4.0  2.9  . 2655316E+06 

SOLUTION  VECTOR 

•  . 1356493 E+ 00 

. 5953607E+00 
. 4599146E+06 


TABLE  XXIV 


lUVVIVIU'RVRUfL  ■  Ul 


IV  UH 


I# 


RESULTS  OF 
EIGHT 


ZUNI  DATA  SEARCH  USING 
EQUAL  MASS  GROUPS 


i 


i 


MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

3.0  3.1  . 1786807E+06 

SOLUTION  VECTOR 
. 1293  27  5E+00 
. 6663778 E+ 00 
. 3094845E+06 


*  MEDIAN  RADIUS  (um)  BETA  STANDARD  DEVIATION 

2.0  3.3  . 18  56012E+06 

SOLUTION  VECTOR 
. 1151535E+00 
. 7367282E+00 


< 


TABLE  XXIV — Continued 


+ 

MEDIAN  RADIUS 
0.6 

(um) 

BETA 

3.9 

STANDARD  DEVIATION 
. 158  516  3E+06 

SOLUTION  VECTOR 
. 9413213E-01 
.  534  5019E+00 
.  2745585E+06 

• 

MEDIAN  RADIUS 
0.5 

(um) 

BETA 

4.0 

STANDARD  DEVIATION 
.  2620418E+06 

• 

SOLUTION  VECTOR 
. 9349972E-01 
. 5644583E+00 
. 4538701E+06 

t 


* 
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sensitivity  analysis,  the  total  mass  of  material  lofted 
by  the  weapon  has  no  effect  on  the  optimum  particle  size 
distribution  that  the  SEARCH4  program  selects.  However,  it 
does  affect  each  of  the  solution  vector  elements  by  the 
same  amount.  Therefore,  in  order  to  make  the  ZUNI  output 
solution  vectors  directly  comparable  to  the  solution  vec¬ 
tors  in  the  model  sensitivity  study,  rather  than  increasing 
each  sample  mass  by  the  same  amount,  the  total  mass  lofted 
was  decreased  by  five  orders  of  magnitude.  This  change 
does  produce  solution  vectors  that  are  directly  comparable 
to  the  ones  in  the  sensitivity  analysis  section  of  Chapter 
III.  Moreover,  it  does  not  affect  the  optimum  particle 
size  distribution  selected  by  the  SEARCH4  program. 

Analysis  of  ZUNI  Airborne 
Sample  Search  Results 

Initially,  the  results  of  the  ZUNI  search  contained 
in  Tables  XIX  through  XXIV  appear  to  be  inconclusive.  In 
all  cases,  for  the  optimum  particle  size  distribution 
selected  by  the  SEARCH4  program,  the  accompanying  solution 
vector  has  two  elements  that  are  within  one  order  of  magni¬ 
tude  of  each  other  and  a  third  element  that  is  six  orders 
of  magnitude  larger  than  the  first  two.  An  acceptable 
solution  vector  would  be  one  with  all  of  its  elements 
equal,  or  at  least,  same  order  of  magnitude.  Without  any 
further  knowledge  of  the  program's  sensitivity  to  certain 
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input,  the  output  contained  in  Tables  XIX  through  XXIV 
does  not  indicate  an  acceptable  solution. 

However,  in  light  of  the  SEARCH4  program  sensi¬ 
tivity  analysis  conducted  in  Chapter  III,  it  is  possible 
that  the  results  of  the  ZUNI  search  can  be  interpreted. 

It  was  discovered  during  the  sample  mass  portion  of  the 
sensitivity  analysis  that  if  two  of  the  sample  masses  were 
in  error  by  as  little  as  10  percent,  the  resultant  program 
output  solution  vector  would  contain  two  elements  that  were 
the  same  order  of  magnitude,  and  one  element  that  was  about 
six  orders  of  magnitude  larger  than  the  other  two.  Yet, 
even  though  the  solution  vector  was  unacceptable,  the  log- 
slope  of  the  optimum  particle  size  distribution,  selected 
by  the  SEARCH 4  program  was  surprisingly  close  to,  if  not 
identical  to  the  actual  log-slope  used  to  construct  the 
input  sample  mass  data. 

Since  little  is  known  about  the  sampling  conditions 
for  this  ZUNI  sample  data  other  than  sample  altitude  and 
time  of  collection,  it  is  most  likely  that  these  samples 
were  not  taken  in  accordance  with  the  standard  sampling 
conditions  assumed  by  the  SEARCH4  program.  For  example, 
the  sampling  aircraft  probably  did  not  fly  through  the  cloud 
center  at  each  sampling  altitude.  Additionally,  each  total 
sample  collection  time  is  longer  than  the  time  that  it 
would  have  taken  the  aircraft  to  fly  through  the  cloud 
once.  This  could  indicate  multiple  passes  at  each  altitude. 


These  differences  could  easily  account  for  a  10  percent 
variation  in  the  sample  masses. 

The  resultant  hypothesis  is  that  because  of  non¬ 
standard  sampling  methods,  two  of  these  ZUNI  airborne 
samples  are  10  percent  greater  than  what  they  should  have 
been.  If  this  is  true,  then  the  lesson  learned  in  the  pro¬ 
gram  sensitivity  analysis  in  Chapter  III  can  be  applied. 
That  is,  even  though  the  solution  vector  is  unacceptable, 
the  optimum  particle  size  distribution  log-slope  selected 
by  the  SEARCH4  program  may  indicate  the  true  log-slope  of 
the  stabilized  nuclear  cloud  particle  size  distribution. 

In  this  case,  the  optimum  log-slope  selected  by 
the  SEARCH 4  program  varies  from  2.7  to  5.0  depending  on 
how  many  equal  mass  groups  are  used  for  the  analysis.  For 
the  ZUNI  analysis,  dividing  the  mass  aloft  into  eight  equal 
mass  groups  and  then  using  only  the  first  three  groups 
assumes  that  63  percent  ((1/8) *5)  of  the  mass  lofted  by 
the  weapon  has  fallen  below  the  sampling  altitudes  during 
the  three  hours  prior  to  sampling.  This  assumption  is 
probably  not  completely  valid,  but  serves  as  a  limiting 
case.  Regardless  of  how  many  groups  were  used  for  the 
analysis,  the  SEARCH4  program  most  frequently  selected 
log-slopes  that  were  between  2.9  and  3.9.  However,  this 
is  by  no  means  a  clear  indication  of  the  actual  particle 
size  distribution  log-slope.  Perhaps  all  that  can  be  said 
at  this  point  is  that  the  integral  method  of  airborne 


V.  The  Differential  Method 


Background 

One  benefit  of  using  the  integral  method  of  air¬ 
borne  filter  sample  analysis  is  that,  aside  from  the 
sample  altitude,  time  of  collection,  and  aircraft  flight 
path  information,  the  total  mass  on  each  sample  is  the  only 
other  piece  of  data  that  is  required.  No  information  con¬ 
cerning  the  particle  sizes  found  in  the  samples  is  needed. 
Based  on  Colarco's  (6:67)  particle  fall  time  versus  radius 
data,  good  assumptions  concerning  the  particle  sizes  that 
should  be  found  in  the  samples  can  be  made.  However,  as 
demonstrated  in  Chapter  IV,  the  integral  method  solutions 
may  not  be  conclusive.  Consequently,  the  theory  for  an 
alternative  method,  the  differential  method,  has  been 
developed. 

This  differential  method  uses  the  relative  number 
of  particles  of  each  size  found  in  the  airborne  sample  in 
conjunction  with  the  sample  altitude,  time  of  collection, 
and  aircraft  flight  path  information  to  unfold  the  stabi¬ 
lized  nuclear  cloud  particle  size  distribution.  Reduced 
airborne  sample  data  from  actual  nuclear  clouds  containing 
the  relative  number  of  particles  of  each  size  on  each  sam¬ 
ple  was  not  available  when  this  study  was  completed.  Hope¬ 
fully,  such  data  will  be  available  for  a  future  study.  In 


anticipation  of  this,  the  differential  method  theory  is 
presented  in  this  chapter. 

Differential  Method  Theory 

Essentially,  the  differential  method  uses  the  same 
particle  fall  mechanics  and  initial  stabilized  cloud  model 
as  presented  in  Chapter  II.  The  major  difference  is  that 
an  assumed  particle  mass-size  distribution  is  not  used. 
Rather,  an  assumed  particle  number-size  distribution  is 
used  to  calculate  the  number  of  particles,  of  each  size 
found  in  the  airborne  sample,  that  would  have  been  present 
at  cloud  stabilization  had  the  assumed  particle  size  dis¬ 
tribution  represented  the  nuclear  cloud  at  stabilization. 
Then,  each  of  these  groups  is  represented  by  a  vertical 
gaussian  and  allowed  to  fall  from  its  initial  lofted  alti¬ 
tude  for  a  time  equal  to  the  sample  collection  time. 
Finally,  the  number  of  particles  in  each  size  group  at  the 
sample  altitude  is  compared  to  the  number  of  particles  of 
that  same  size  found  on  the  sample.  In  other  words,  a 
ratio  is  made  of  the  calculated  number  of  particles  of  a 
given  size  at  a  given  altitude  to  the  measured  number  of 
particles  in  an  airborne  sample  of  that  same  size  at  that 
same  altitude.  If  all  of  these  ratios  for  all  the  groups 
are  equal  or  close  to  equal,  then  the  assumed  particle 
size  distribution  is  the  one  that  must  have  represented 
the  nuclear  cloud  at  stabilization.  One  obvious  advantage 
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to  this  method  is  that  only  one  good  airborne  sample  is 
required  for  the  analysis. 

The  following  algorithm  outlines  the  differential 
method  of  airborne  sample  analysis: 

1.  Airborne  sample  data,  to  include  sample  alti¬ 
tude,  time  )f  collection,  and  relative  number  of  particles 
per  particle  radius  found  on  the  sample  are  input. 

2.  The  first  assumed  lognormal  particle  size  dis¬ 
tribution  is  used  to  calculate  the  number  of  particles  in 
each  of  the  measured  sizes. 

3.  A  vertical  gaussian  representing  each  of  these 
groups  is  allowed  to  fall  from  its  initial  lofted  altitude 
for  a  time  equal  to  the  sample  collection  time. 

4.  The  number  of  particles  at  the  sample  altitude 
from  each  of  these  falling  groups  is  calculated. 

5.  This  calculated  number  of  particles  is  com¬ 
pared  to  the  measured  number  of  particles  in  the  airborne 
sample  in  ratio. 

6.  All  the  ratios  from  all  the  different  size 
groups  are  compared  to  each  other  and  a  standard  deviation 
from  an  algebraic  mean  for  these  ratios  is  calculated. 

7.  Another  trial  distribution  is  selected  and  the 
process  is  begun  again.  The  trial  distribution  that  pro¬ 
duces  the  smallest  standard  deviation  among  the  group 
ratios  is  the  particle  size  distribution  that  best  repre¬ 
sents  the  stabilized  nuclear  cloud. 


This  differential  method  of  airborne  sample  analy¬ 
sis  may  seem  simpler  than  the  integral  method  since  it 
requires  only  one  sample.  However,  there  are  some  limita¬ 
tions  to  this  method.  For  example,  there  is  considerable 
difficulty  involved  in  obtaining  an  accurate  breakdown  of 
the  relative  number  of  particles  of  each  size  on  an  air¬ 
borne  sample  (15:21).  Additionally,  the  same  corrections 
for  nonstandard  sampling  conditions  that  applied  to  the 
integral  method  apply  to  the  differential  method  as  well. 
Moreover,  rather  than  base  conclusions  about  the  nature  of 
a  nuclear  cloud  particle  size  distribution  on  the  results 
of  one  sample  analysis,  it  is  recommended  that  the  integral 
method  of  airborne  sample  analysis  be  used  in  conjunction 
with  the  differential  method.  Then  the  results  from  both 
methods  of  analysis  can  be  used  to  predict  the  stabilized 
nuclear  cloud  particle  size  distribution. 
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VI .  Conclusions  and  Recommendations 

Conclusions 

In  order  to  accurately  predict  fallout  patterns, 
the  nuclear  cloud  particle  size  distribution  must  be  known. 
However,  analysis  of  fallen  material  alone  may  lead  to 
inaccurate  conclusions  concerning  the  particle  size  dis¬ 
tribution  in  the  nuclear  cloud.  For  example,  in  Hopkins' 
analysis  of  the  Mt.  St.  Helens  ash  cloud  he  discovered  that 
because  of  breakage  on  impact  and  sieving  of  the  fallen 
material,  the  particle  size  distribution  calculated  from 
the  fallen  material  analysis  could  not  be  used  to  model 
the  observed  fallout  pattern  (17:86).  Likewise,  for 
nuclear  clouds,  analysis  of  fallen  material  alone  may  be 
misleading.  Therefore,  the  best  estimate  of  the  actual 
falling  particle  size  distribution  in  a  nuclear  cloud  can 
be  made  from  analysis  of  fallen  material  samples  as  well 
as  analysis  of  cloud  samples  taken  at  different  altitudes 
and  different  times.  This  study  presented  two  methods  for 
unfolding  the  nuclear  cloud  particle  size  distribution  from 
nuclear  cloud  sample  analysis. 

The  integral  method  of  nuclear  cloud  sample  analysis 
presented  in  Chapter  II  is  a  viable  method  of  unfolding  the 
stabilized  nuclear  cloud  particle  distribution  if  the 
sampling  conditions  were  well  documented.  If  not,  this 


method  produces  results  that  are  inconclusive.  The  major 
advantage  of  the  integral  method  is  that,  aside  from  the 
standard  information  concerning  where  and  when  the  samples 
were  taken  with  respect  to  the  nuclear  cloud  location,  the 
only  other  sample  data  required  is  the  total  mass  on  each 
of  the  filter  samples.  Additionally,  if  total  activity  on 
each  of  the  filter  samples  at  a  given  time  after  the  burst 
is  known,  the  integral  method  of  activity  summation  rather 
than  mass  summation  can  be  used  to  unfold  the  stabilized 
nuclear  cloud  particle  distribution  if  some  additional 
assumptions  are  made.  These  assumptions  are  incorporated 
into  the  code  by  minor  modifications  that  include  using 
Freiling's  2.5  moment  approximation  (12:6)  for  the  activity 
distribution  versus  particle  size,  and  the  Way-Wigner 
(23:1318)  approximation  for  fission  product  decay. 

If  a  more  detailed  reduction  of  sample  data, 
including  the  relative  number  of  particles  of  each  size 
on  each  of  the  filter  samples  is  available,  the  differen¬ 
tial  method  of  nuclear  cloud  sample  analysis  presented  in 
Chapter  IV  can  be  used  to  unfold  the  nuclear  cloud  par¬ 
ticle  size  distribution.  Unlike  the  integral  method,  the 
differential  method  does  not  require  as  much  information 
concerning  the  actual  sampling  conditions,  aside  from 
sample  altitude  and  sample  time  of  collection.  For  example, 
since  particles  of  all  sizes  in  a  nuclear  cloud  are  assumed 
to  be  normally  distributed  in  the  horizontal  plane,  the 


details  on  the  cloud  penetration  other  than  altitude  of 
penetration  are  not  important.  Only  the  relative  number 
of  particles  per  radius  is  important  and  this  will  not  vary 
with  different  aircraft  penetration  points  at  any  given 
altitude.  However,  the  tradeoff  is  that  the  more  detailed 
sample  data  reduction  for  the  early  U.S.  atmospheric  tests 
is  difficult  to  acquire. 

The  optimum  approach  to  unfolding  the  nuclear  cloud 
particle  size  distribution  from  airborne  cloud  samples  is 
to  acquire  a  set  of  filter  sample  data  that  will  allow  the 
use  of  both  the  integral  and  differential  methods  presented 
in  this  study.  Obviously,  if  enough  data  is  available  for 
a  given  nuclear  cloud  to  allow  the  use  of  both  of  these 
methods,  then  the  cloud  particle  size  distribution  can  be 
predicted  with  greater  certainty. 

Recommendations 

There  are  three  recommendations  to  be  made.  The 
first  one  is  that  a  thorough  search  be  made  of  the  classi¬ 
fied  archives  at  each  of  the  national  laboratories,  Lawrence 
Livermore  and  Los  Alamos,  for  reduced  cloud  sample  data. 

Both  labs  were  sponsors  for  the  atmospheric  tests  and  there 
is  a  high  probability  that  sufficiently  reduced  sample  data 
already  exists  in  their  archives.  In  most  cases,  the  type 
of  data  required  by  the  two  methods  presented  in  this  study 
is  unclassified.  Such  was  the  case  for  the  ZUNI  data  used 
in  Chapter  IV. 


The  remaining  two  recommendations  concern  modifica¬ 
tions  to  the  code  contained  in  Appendix  B.  A  better  method 
should  be  developed  for  conducting  the  search  rather  than 
trial  and  error.  Also,  the  code  should  be  modified  to 
accept  the  required  data  and  perform  both  integral  and  dif¬ 
ferential  analysis  of  airborne  samples.  These  modifica¬ 
tions,  along  with  a  set  of  sampling  standards  could  be  used 
to  facilitate  the  analysis  of  data  from  any  future  sampling 
endeavors . 


Appendix  A:  Empirical  Cloud  Loft  Model 


This  appendix  contains  the  empirical  equations 
used  to  spatially  position  the  stabilized  nuclear  cloud. 

Two  sets  of  equations  are  used.  The  first  set  was  devel¬ 
oped  by  Hopkins  using  a  polynomial  least  squares  fit  from 
DELFIC  output  for  various  yields.  With  weapon  yield  in 
kilotons  and  particle  radius  in  micrometers  as  input, 
these  equations  can  be  used  to  calculate  the  altitude  of 
a  vertical  gaussian  representing  any  mono-size  group  of 
particles  at  cloud  stabilization  time  (16:14-15): 

H  =  (SLOPE  *  2  *  R)  +  (INTERCEPT) 

9 

SLOPE  =  -  EXP {1.574  -  0.01197* In (YKT)  +  0 . 03636*ln (YKT) 2 
-  0 . 0041*ln (YKT) 3  +  0 . 0001965*ln (YKT) 4} 

INTERCEPT  =  EXP {7.889  +  0.34*ln(YKT)  +  0 . 001226*ln (YKT) 2 
-  0.005227* In (YKT) 3  +  0 . 000417*ln (YKT) 4 } 

where 

H  =  the  average  lofted  altitude,  in  meters,  of 
g  the  center  of  a  vertical  gaussian  representing 
a  mono-size  group  of  particles  of  radius  R 

R  =  the  radius  in  micrometers  of  particles  in  the 
mono-size  group 

YKT  =  the  weapon  yield  in  kilotons 


The  next  set  of  empirical  equations  were  developed 
by  Conners  using  methods  similar  to  those  used  by  Hopkins, 
above.  These  equations  are  used  to  calculate  the  standard 
deviation  for  each  of  the  mono-size  group  vertical  gaussian 
distributions  (7:19-20): 


Zg  =  (SLOPE  *  2  *  R)  +  (INTERCEPT) 


SLOPE  =  7  -  EXP(l. 78999  -  0.048249* In (YKT) 
+  0.0230248* In (YKT) 2  -  0 . 00225965*ln (YKT) 3 


where 


+  0.000161519* In (YKT) 4 } 


INTERCEPT  =  EXP{ 7 . 03518  +  0 . 158914*ln (YKT) 
+  0.0837539* In (YKT) 2  -  0 . 0155464*ln (YKT) 3 
+  0.00086 21 03* In (YKT) 4) 


=  the  predicted  vertical  thickness  in  meters  of 
a  mono-sized  group  of  particles  with  radius  R 

=  the  radius  in  micrometers  of  particles  in  the 
mono-size  group 

=  the  weapon  yield  in  kilotons 


If  Zg  is  assumed  to  be  a  two-sigma  distribution, 
then  the  standard  deviation,  in  meters,  for  each  mono-size 
group  vertical  gaussian  is  defined: 


o  =  0.25  *  Z 
9  9 
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Appendix  B:  SEARCH 4  Computer  Code  and  Glossary 

of  Program  Variables 


This  appendix  contains  the  computer  code  that 


employs  the  integral  method  of  airborne  filter  sample 


analysis.  It  is  written  in  FORTRAN-5  (FORTRAN-77)  and  was 


run  on  a  CYBER  170-845.  Following  the  code  listing  are 


example  input  files,  FILTER  and  SEARCH,  and  a  glossary  of 


all  the  variables  used  in  the  program. 


t**********************  SEARCH4  ********************** 


FORTRAN- 7 7  (FORTRAN- 5)  ON  CYBER 


**  13  NOVEMBER  1985  **  A  FORTRAN  VERSION  OF  FILTER4C  * 


**  (BASIC  PROGRAM)  THAT  SEARCHES  FROM  BETA=2  to  BETA=5  * 


FOR  AN  APPROPRIATE  SOLUTION  VECTOR.  PARTICLE 


**  SIZE  RANGE  ON  THE  FILTER  SAMPLES  IS  FROM  A  FEW 


MICRONS  TO  70-120  MICRONS. 


INPUT:  FROM  FILES  "FILTER"  AND  "SEARCH" 


ALSO,  TWO  INTERACTIVE  QUESTIONS. 


OUTPUT:  TO  A  FILE  "ANSWER" 


**  TYPE  STATEMENTS  ** 


PROGRAM  START 


REAL  INMASS (25),  ZZERO(25),  R(25),  SIGMA(25),  TOTMAS(25) 


REAL  TAME (25),  ALTCHK (25)  ,  A(25,25),  B(25,50),  C(25) 
REAL  RHOF,  YKT,  MASS,  TONE,  BETA,  ALPHAl ,  ISTEP,  G 
REAL  ALPHA2,  NUMGRP,  MOMENT,  RAT,  BP,  PD INC,  PDCHK 
REAL  Z,  PD,  RIGHT,  LEFT,  HOLT,  RHO,  AIDA,  TEMPAL,  R2CD 
REAL  REY,  VZ ,  SLPDRG,  DELTAT,  ALTZ ,  DZ ,  DUMMY,  SUMMER 
REAL  DTHOLD ,  GAUSS,  BF,  TEMP,  TOTAL,  DENOM 
REAL  YOU,  WIN,  CRAB,  AVERAG,  FISH,  PIG 
REAL  EXTRAS 

INTEGER  PLUTO,  I,  MARK,  COLT,  FLAG,  J,  T,  FIGS,  OMEGA 
INTEGER  M,  K,  N,  Q,  P,  H,  F,  U,  X,  Y,  MINUSE,  QUICK 
INTEGER  ZETA 
INITIALIZATION  ** 

RHOF  =  2600.0 
G  =  9.77 
MARK  =  0 
RIGHT  =  0 
LEFT  =  0 
FLAG  =  0 
PDCHK  =  0 
COLT  =  0 
FIGS  =  0 
HOLT  =  0 
SUMMER  =  0 
CRAB  =  0 
MINUSE  =  0 


FISH 


0 


OMEGA  =  0 


ZETA  =  0 

**  READ  FILTER  DATA  INPUT  FILE  ** 

OPEN (UNIT=7,FILE=' FILTER' ) 

REWIND  7 
READ (7, 107)  YKT 
PRINT*,  'YKT  =  ' , YKT 
READ (7 ,111)  PLUTO 

PRINT*,  'NUMBER  OF  SAMPLES  =  ', PLUTO 
111  FORMAT (16) 

READ (7, 107)  MOMENT 

PRINT*,  'DISTRIBUTION  MOMENT  =  ', MOMENT 

107  FORMAT (F6.0) 

READ (7, 108)  MASS 

PRINT*,  'MASS  ALOFT  =  ' ,MASS 

108  FORMAT (E13. 2) 

DO  10  I  =  1,  PLUTO 

READ ( 7 , 109 )  TOTMAS ( I ) , ALTCHK ( I ) , TONE 
PRINT*,  'SAMPLE  NUMBER  ',1,'  :  MASS  =  ' 

$ ALTITUDE  =  ', ALTCHK (I) , '  TIME  =  * , TONE 

109  FORMAT (F9.0,F5.0,F5.1) 

TAME  (I)  =  TONE*60 . 0 

10  CONTINUE 

OPEN (UNIT=6,FILE=' SEARCH ' ) 


TOTMAS ( I ) , ' 


REWIND  6 

READ (6 ,349)  BETA, ALPHA1 , ALPHA2 , ISTEP ,NUMGRP 


349  FORMAT (F4 .1,F6 .1,F7. 2, F 5.1 ,F5.1) 
**  INTEROGATE  HUMAN  FOR  INPUT  ** 


PRINT*,  '  ' 

PRINT*,  'DO  YOU  WANT  A  QUICK  SEARCH  ?' 

PRINT*,  ' (ENTER  1  FOR  QUICK  SEARCH) ' 

READ*,  QUICK 
IF  (QUICK  .EQ.  1)  THEN 
DENOM  =  1.0 

ELSE 

DENOM  =10.0 
END  IF 

I STEP  =  I STEP /DENOM 
PRINT*,  '  ' 

PRINT*,  'DO  YOU  WANT  TO  USE  MORE  EQUAL  MASS  GROUPS  ?' 

PRINT*,  'IF  SO,  ENTER  THE  NUMBER  TO  BE  ADDED  TO  NUMGRP.' 

PRINT*,  'IF  NOT,  ENTER  0.’ 

READ*,  EXTRAS 

NUMGRP  =  NUMGRP  +  EXTRAS 

PRINT*,  '  ' 

PRINT*,  'ENTER  THE  OUTPUT  PARAMETER.' 

PRINT*,  ' (ENTER  1  FOR  SOLUTION  VECTORS  WITH  MINUSES  <=  1) ' 

PRINT*,  '(ENTER  0  FOR  SOLUTION  VECTORS  WITH  MINUSES  =  0)' 

READ*,  ZETA 
IF  (ZETA  .NE.  1)  THEN 
ZETA  =  0 
END  IF 
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**  END  OF  INPUT  ** 

PRINT*,'  CURRENTLY  WORKING  ON  BETA  =  ' ,BETA 

**  OPEN  OUTPUT  FILE  AND  WRITE  HEADINGS  ** 

OPEN (8 ,FILE= 'ANSWER' , STATUS= 'NEW ' ) 

101  FORMAT ('  ') 

**  GOSUB  FOR  PARTICLE  DISTRIBUTION  CALCULATION  ** 

1  CONTINUE 

DO  15  I  =  1,  PLUTO 

CALL  PARTY (R(I) , BETA , ALPHAl , NUMGRP , MOMENT , I STEP , 

$ ALPHA 2 , COLT , PLUTO , FLAG , PDCHK , OMEGA , DENOM , EXTRAS ) 

IF  (OMEGA  .EQ.  1)  THEN 
GO  TO  35 
END  IF 

15  CONTINUE 

DO  20  I  =  1,  PLUTO 
INMASS (I)  =  MASS /NUMGRP 
20  CONTINUE 

**  GOSUB  FOR  GROUP  LOFTED  ALTITUDE  AND  SIGMA  CALCULATION  ** 

DO  30  I  =  1, PLUTO 

CALL  HOPKIN (R ( I ) ,ZZERO(I) , SIGMA (I) ,YKT, PLUTO) 

30  CONTINUE 

**  GOSUB  TO  USAIR  WITH  THE  LARGEST  PARTICLE'S  LOFTED  ALTITUDE  ** 
TEMPAL  =  ZZERO (PLUTO) -700.0 
CALL  USAIR (TEMPAL, RHO, AIDA) 

**  CALCULATE  DELTAT  ** 

R2CD  =  (32. 0  *RHO*RHOF*G* (R (PLUTO) **3 ) ) / (3.0* (AIDA**2) ) 
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IF  (R2CD  .LT.  120.0)  THEN 


* 


& 


REY  =  (R2CD/24 .0) - (2.3363E-4* (R2CD**2) )+ (2 .0154E 
$-6* (R2CD**3) ) - ( 6 . 9105E-9* (R2CD**4 ) ) 

ELSE 

REY  =  10.0**((-1.29536)+(.986*ALOG10(R2CD))- 
$ ( .046677* (ALOGIO (R2CD) ) **2 ) + ( 1 . 1235E-3* (ALOGIO (R2CD) ) 
$**3 )  ) 

END  IF 

VZ  =  (REY*AIDA) / ( 2 . 0*RHO*R (PLUTO) ) 

SLPDRG  =  1 . 0+ {1 . 16  5E-7/ (R (PLUTO) *RHO) ) 

DELTAT  =  1400 .0/ (VZ*SLPDRG) 

**  MAIN  PROGRAM======>>  PARTICLE  GROUP  FALL  MECHANICS  LOOPS 

$  <<===  ** 

DO  50  J  =  1, PLUTO 
DO  60  T  =  1, PLUTO 
ALTZ  =  ZZERO(T) 

36  CONTINUE 

CALL  USAIR ( ALTZ, RHO, AIDA) 

R2CD  =  (32 . 0*RHO*RHOF*G* (R(T)**3))/(3.0* (AIDA** 2) ) 

IF  (R2CD  .LT.  120.0)  THEN 

REY  =  (R2CD/ 24.0)-(2.3363E-4* (R2CD**2 ) ) + ( 2 . 0154E 
$-6* (R2CD**3 ) ) - (6 . 9105E-9* (R2CD**4) ) 

ELSE 

REY  =  10.0**( (-1.29536)+(. 986*ALOG10 (R2CD) ) - 
$ ( .046677* (ALOGIO (R2CD) ) **2) + (1 . 1235E-3* (ALOGIO ( R2CD) ) 


4. 

K 


END  IF 

IF  (ALTZ  .LT.  0)  THEN 
GO  TO  31 
END  IF 

VZ  =  (REY*AIDA) / ( 2 . 0*RHO*R (T) ) 

SLPDRG  =  1 . 0+ (1 . 165E-7/ (R  (T) *RHO) ) 

DZ  =  VZ* SLPDRG *DELTAT 
31  CONTINUE 

ALTZ  =  ALTZ  -  DZ 
DUMMY  =  INMASS (T) 

IF  (ALTZ  .LE.  0)  THEN 
GAUSS  =  0 
GO  TO  32 
END  IF 

GAUSS  =  (DUMMY*EXP (-. 5* (( (ALTCHK(J) -ALTZ) /SIGMA (T) ) 
$  **  2 ) ) ) / (SQRT (2.0*3.14159) * SIGMA (T) ) 

SUMMER  =  SUMMER  +  DELTAT 
IF  (GAUSS  .LT.  1.0E-10)  THEN 
GAUSS  =  0 
GO  TO  32 
END  IF 

IF  (FIGS  .EQ.  1)  THEN 
GO  TO  34 
END  IF 

IF  ( (SUMMER+DELTAT)  .GT.  TAME (J) )  THEN 
FIGS  =  1 
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DTHOLD  =  DELTAT 


DELTAT  =  TAME ( J )  -  SUMMER 
GO  TO  36 

ELSE 

GO  TO  36 
END  IF 

34  CONTINUE 

DELTAT  =  DTHOLD 
FIGS  =  0 
32  CONTINUE 

A ( J , T)  =  GAUSS 
SUMMER  =  0 

60  CONTINUE 
50  CONTINUE 

**  FINAL  DETERMINATION  OF  THE  A  VALUES  AND  A  VALUE  COMPARISON  ** 

**  1.  INVERSION  OF  THE  COEFFICIENT  MATRIX  ** 

DO  61  T  =  1, PLUTO 
DO  62  H  =  1, PLUTO 
A (T, H)  =  A ( T , H ) *NUMGRP 
62  CONTINUE 

61  CONTINUE 
N  =  PLUTO 

**  SET-UP  OF  WORKAREA  MATRIX  WITH  IDENTITY  MATRIX  ON  THE  RIGHT  ** 
DO  63  I  =  1,N 
DO  64  J  =  1 , N 
B  ( I , N+J)  =  0 
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B  ( I ,  J)  =  A  ( I ,  J) 


64  CONTINUE 

B (I , I+N)  =  1 
63  CONTINUE 

**  GAUSS-JORDIN  ELIMINATION  TO  DETERMINE  THE  INVERSE 
DO  65  K  =  1 , N 
IF  (K  .EQ.  N)  THEN 
GO  TO  53 
END  IF 
M  =  K 

DO  67  I  =  (K+l ) , N 
YOU  =  ABS (B ( I / K) ) 

WIN  =  ABS (B (M,K) ) 

IF  (YOU  .GT.  WIN)  THEN 


M  =  I 


END  IF 


67  CONTINUE 


IF  (M  .EQ.  K)  THEN 


GO  TO  53 


END  IF 


DO  68  J  =  K, ( 2*N) 
BF  =  B  (K,  J) 

B  (K,  J)  =  B  (M,  J) 

B (M , J )  =  BF 


68  CONTINUE 


53  CONTINUE 


IF  (B(K,K)  .EQ.  0)  THEN 


I 


# 


F 


GO  TO 

89 

END 

IF 

B  (K 

,J)  =  B 

(K,J)/B(K,K) 

66 

CONTINUE 

IF 

(K  .EQ. 

1)  THEN 

GO  TO 

69 

END 

IF 

DO 

70  I  = 

1, ( K— 1 ) 

DO 

71  J  = 

(K+l) , ( 2*N) 

B(I 

,  J)  =  B 

(I,J)  -  B ( I , 

71 

CONTINUE 

70 

CONTINUE 

IF 

(K  .EQ. 

N)  THEN 

GO  TO  72 


69 

END  IF 

CONTINUE 

DO  73  I  = 

(K+l)  ,N 

DO  74  J  = 

(K+l) , (2*N) 

B(I,J)  = 

B  ( I ,  J )  -  B  ( I 

74 

CONTINUE 

73 

CONTINUE 

65 

CONTINUE 

72 

CONTINUE 

K)  *B  (K,  J) 


K)  *B (K, J) 


**  RETRIEVE  THE  INVERSE  FROM  THE  RIGHT  SIDE  OF  B 


ft 

* 

i 


DO  75  I  =  1 ,N 
DO  76  J  =  1 , N 
B ( I , J)  =  B (I ,  J+N) 

76  CONTINUE 
75  CONTINUE 

**  SHIFT  B  TO  A  FOR  NEXT  OPERATION  ** 

DO  77  Q  =  1 , N 
DO  78  P  =  1 , N 
A  (Q,  P)  =  B  (Q ,  P ) 

78  CONTINUE 

77  CONTINUE 

**  2.  MULTIPLY  A,  THE  INVERSE  MATRIX  TIMES  B,  THE  SAMPLE  VECTOR  ** 
M  =  N 
I  =  N 
J  =  1 
TOTAL  =  0 

**  INPUT  THE  RIGHT  MATRIX  ** 


DO  79 

F  = 

1,1 

DO  80 

U  = 

1,J 

B  (F,  U) 

=  TOTMAS 

(F) 

80 

CONTINUE 

79 

CONTINUE 

DO  81 

I  = 

1  ,M 

DO  8  2 

X  = 

1,  J 

DO  83 

Y  = 

1,N 

TEMP  = 

:  A  ( I 

, Y)  *B 

(Y,  X) 

TOTAL  =  TOTAL  +  TEMP 


IF  (Y  .EQ.  N)  THEN 
C (I)  =  TOTAL 
TOTAL  =  0 
END  IF 

83  CONTINUE 
82  CONTINUE 
81  CONTINUE 

**  3.  FINALLY,  A  LOGNORMAL  CURVE-FIT  CHECK  CODE  ** 
DO  87  I  =  1, PLUTO 
CRAB  =  CRAB  +  C (I) 

IF  (C  ( I )  .LT.  0)  THEN 

MINUSE  =  MINUSE  +1 
IF  (MINUSE  .GT.  ZETA)  THEN 
GO  TO  89 
END  IF 
END  IF 

87  CONTINUE 

PIG  =  PLUTO 
AVERAG  =  CRAB /PIG 
DO  91  I  =  1, PLUTO 
FISH  =  FISH  +  (C (I) -AVERAG) **2 
91  CONTINUE 

FISH  =  SQRT (FISH/ (PIG-1) ) 

IF  (FISH  .GT.  1000000)  THEN 
GO  TO  89 

98 


* 


# 


END  IF 

HOLT  =  ALPHA1-ISTEP 
WRITE (8,101) 

WRITE (8, 101) 

/ 

WRITE  (8,100) 

100  FORMAT ('  ALPHA  BETA  STANDARD  DEVIATION 

$  MINUSES ' ) 

WRITE(8 ,228)  HOLT, BETA, FISH ,MINUSE 
228  FORMAT (2X,F8.5,4X,F5.1,7X,E15.7,20X,I2) 

WRITE (8,101) 

WRITE (8, 555) 

555  FORMAT (29X, 'SOLUTION  VECTOR* ) 

DO  93  I  =  1, PLUTO 
WRITE (8,371)  C  (I ) 

371  FORMAT (22X, El8 .7) 

93  CONTINUE 

WRITE (8 , 101) 

89  CONTINUE 
CRAB  =  0 
FISH  =  0 
AVERAG  =  0 
MINUSE  =  0 

**  GO  GET  ANOTHER  DISTRIBUTION  TO  CHECK  ** 

GO  TO  1 
35  CONTINUE 
REWIND  8 
END 

99 


*************************************** 


# 


I 

■ 


# 


**  SUBROUTINE  PARTY  ** 

★****★***★★*★★*★****★*★*★**★★★★★★****★★ 

SUBROUTINE  PARTY (ROOT, BETA, ALPHAl , NUMGRP , MOMENT, ISTEP, 

| 

$ ALPHA2 , COLT , PLUTO  f  FLAG , PDCHK , OMEGA , DENOM , EXTRAS ) 

REAL  ROOT, BETA, ALPHAl, NUMGRP, MOMENT, I STEP, ALPHA2, 

$ALPHA3 , RAT , BP 

REAL  R (25) , PD INC , PDCHK , PD , Z , LEFT, RIGHT, DENOM 
REAL  EXTRAS 

INTEGER  PLUTO, FLAG, COLT, I, OMEGA 
IF  (COLT  .EQ.  1)  THEN 

READ (6,349)  BETA , ALPHAl , ALPHA2 , ISTEP , NUMGRP 
349  FORMAT (F4.1,F6.1,F7.2,F5.1,F5.1) 

NUMGRP  =  NUMGRP  +  EXTRAS 

PRINT*,'  CURRENTLY  WORKING  ON  BETA  =  ' , BETA 

ISTEP  =  ISTEP/DENOM 
COLT  =  0 
END  IF 

IF  (BETA  .EQ.  0)  THEN 
OMEGA  =  1 
RETURN 
END  IF 

RAT  -  10000.0 
LEFT  =  0 
RIGHT  =  0 
BP  =  ALOG (BETA) 
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ALPHA3  =  ALOG ( ALPHAl )  +  (MOMENT*BP**2 ) 
IF  (FLAG  .GT.  0)  THEN 


|;i 

! 

r? 


GOTO  5 


END  IF 


PDINC  =  (100 . 0/NUMGRP) / 100 . 0 
PDCHK  =  PD INC/ 2.0 


CONTINUE 


Z  =  (ALOG (RAT) -ALPHA3 ) /BP 


IF (Z  .GE.  0)  THEN 

PD=1 . 0- . 5/ ( 1. 0+. 196 8 54*Z+. 115194*Z* *2+. 000344 *Z 


$**3+.019527*Z**4) *  *4 


ELSE 


Z  =  ABS(Z) 

PD= .  5/  (1.0+.196854*Z+.115194*Z**2+.000344*Z**3+ 


$.019527*Z**4) **4 


END  IF 


P 

I  * 


IF  (PD  .GT.  PDCHK)  THEN 


RIGHT  =  RAT 


RAT  =  (RIGHT-LEFT) / 2. 0  +  LEFT 


ELSE 


LEFT  =  RAT 


RAT  =  (RIGHT-LEFT)/ 2.0  +  LEFT 


END  IF 


IF  ( (ABS (PD-PDCHK) )  .LT.  .00001)  THEN 


PDCHK  =  PDCHK  +  PDINC 


FLAG  =  FLAG  +  1 


1 


ROOT 


RAT  *  .000001 


# 


RAT  =  10000.0 
IF  (FLAG  .GE.  PLUTO)  THEN 
GOTO  6 

ELSE 

GOTO  7 
END  IF 

ELSE 

GOTO  5 
END  IF 

6  CONTINUE 

ALPHAl  =  ALP HA 1  +  I STEP 
IF  (ALPHAl  .GT.  ALPHA2 )  THEN 
COLT  *  1 
END  IF 
FLAG  =  0 

7  CONTINUE 
RETURN 
END 

***************************************** 

**  SUBROUTINE  HOPKIN  ** 

***************************************** 

SUBROUTINE  HOPKIN (RAD, INALT, SIG, YKT, PLUTO) 

REAL  RAD , INALT , S IG , YKT , IM , SM , ID , SD , DELTAZ , LOGYKT 
INTEGER  PLUTO 
LOGYKT  =  ALOG (YKT) 
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IM  =  EXP (7.889  +  . 34*LOGYKT  +  . 001226 *LOGYKT**2 


$ .005227*LOGYKT**3  +  . 000417*LOGYKT**4 ) 

SM  =  -EXP (1.574  -  . 01197*LOGYKT  +  . 03636*LOGYKT**2  - 
$ .0041*LOGYKT**3  +  . 000196 5*LOGYKT**4 ) 

ID  =  EXP (7.03518  +  . 158914*L0GYKT  +  .  0837539*LOGYKT 
$**2  -  . 0155464*LOGYKT**3  +  .  000862103*LOGYKT**4 ) 

SD  =  7.0  -  EXP(1. 78999  -  . 048249*LOGYKT  +  .0230248* 
$LOGYKT**2  -  . 0022596 5*LOGYKT**3  +  .000161519*LOGYKT**4) 
INALT  =  IM  +  (2 . 0*RAD* 1000 000 . 0*SM) 

DELTAZ  =  ID  +  (2 . 0*RAD*1000000 . 0*SD) 

SIG  =  . 25*DELTAZ 

RETURN 

END 

ar***************************************** 

**  SUBROUTINE  USAIR  ** 

SUBROUTINE  USAIR (Z ,RHO, AIDA) 

REAL  Z , RHO, AIDA, LK , TZ , PZ 
IF  (Z  .LT.  11000)  THEN 
LK  =  -.006545 
TZ  =  288.15  -  ( . 006 545 *Z ) 

PZ  =  101300.0  *  (288. 15/TZ) ** (-.134164/. 006545) 


GO  TO  22 


1 


END  IF 


IF  (Z  .LT.  20000)  THEN 


LK  =  0 


XX  M  W* 


m 


> 

5- 


i 


# 


# 


4 


TZ  =  216.65 

PZ  =  22690 . 0*EXP (-.034164* (Z-11000.0) / 216 . 65) 
GO  TO  22 
END  IF 

IF  (Z  .LT.  32000)  THEN 
LK  =  .001 

TZ  =  216.65  +  .001* (Z-20000.0) 

PZ  =  5528*(216.65/TZ)**(. 034164/. 001) 

GO  TO  22 
END  IF 

IF  (Z  .LT.  47000)  THEN 
LK  =  .0028 

TZ  =  228.65+.0028MZ-32000.0) 

PZ  =  888. 8*(228.65/TZ) **(.034164/. 0028) 

GO  TO  22 


•  4 


END  IF 

IF  (Z  .LT.  53000)  THEN 
LK  =  0 
TZ  =  270.65 

PZ  =  115. 8*EXP(-. 034164* (Z-47000.0) /270. 65) 
GO  TO  22 
END  IF 

IF  (Z  .LT.  59000)  THEN 
LK  =  -.0028 

TZ  =  265.0 5-. 0028* (Z-53000.0) 

PZ  =  54. 87* (265.05/TZ) ** (-. 034164/. 0028) 
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GO  TO  22 


END  IF 


IF  (Z  .LT.  70000)  THEN 


LK  =  -.0028 


TZ  =  248. 2 5-. 0028* (Z-59000.0) 


PZ  =  25. 132* (248 .25/TZ) ** (-. 034164/ .0028) 


GO  TO  22 


END  IF 


IF  (Z  .GT.  70000)  THEN 


PRINT*, ’ 


WARNING. . .ALTITUDE  ERROR! ! ! ' 


END  IF 


CONTINUE 


RHO  =  .003484* (PZ/TZ) 


AIDA  =  (1.458E-6* (TZ**1.5) ) / (TZ+110.4) 


RETURN 


EXAMPLE  INPUT  FILES 


FILTER 

3  500  .  YIELD  IN  KILOTONS 

3  .  NUMBER  OF  SAMPLES 

3.0  .  DISTRIBUTION  MOMENT 

1.06E+9  .  MASS  LOFTED  IN  KILOGRAMS 

.255  13106182.0  _  1 

.3118  15758151.5  _ !  SAMPLE:  MASS  (kg),  ALTITUDE  (m) , 

.142  16246184.0  _ ]  AND  TIME  OF  COLLECTION  (min) 


SEARCH 

LOG-SLOPE  STARTING  FINAL  INCREMENT  NUMBER  OF  GROUPS 


MEDIAN  MEDIAN 
RADIUS  RADIUS 


1.5 

100.0 

300.00 

5.0 

3.0 

1.6 

100.0 

300.00 

5.0 

3.0 

1.7 

100.0 

200.00 

5.0 

3.0 

1.8 

100.0 

200.00 

5.0 

3.0 

1.9 

100.0 

200.00 

5.0 

3.0 

2.0 


20.0  100.00 


2.0 


3.0 


2.1 

20.0 

75.00 

1.0 

3.0 

2.2 

20.0 

75.00 

1.0 

3.0 

2.3 

10.0 

40.00 

1.0 

3.0 

2.4 

1.0 

30.00 

1.0 

3.0 

2.5 

1.0 

30.00 

1.0 

3.0 

2.6 

1.0 

30.00 

1.0 

3.0 

2.7 

1.0 

30.00 

1.0 

3.0 

2.8 

1.0 

20.00 

1.0 

3.0 

2.9 

1.0 

20.00 

1.0 

3.0 

3.0 

1.0 

20.00 

1.0 

3.0 

3.1 

1.0 

20.00 

1.0 

3.0 

3.2 

1.0 

20.00 

1.0 

3.0 

3.3 

1.0 

10.00 

0.5 

3.0 

3.4 

1.0 

10.00 

0.5 

3.0 

3.5 

0.5 

6.00 

0.5 

3.0 

3.6 

0.1 

3.00 

0.2 

3.0 

3.7 

0.1 

2.00 

0.2 

3.0 

3.8 

0.1 

2.00 

0.2 

3.0 

3.9 

0.1 

1.00 

0.1 

3.0 

4.0 

0.1 

1.00 

0.1 

3.0 

4.1 

0.1 

0.60 

0.1 

3.0 

4.2 

0.1 

0.40 

0.1 

3.0 

4.3 

.110 

.200 

.180 

3.0 

4.4 

.090 

.150 

.150 

3.0 

4.5 

.080 

.120 

.120 

3.0 
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KAWIf  1  I' 


Variable 

Name  Real/ Integer 

INMASS (25)  Real 
Z  Z  ERO (25) 

R  ( 25 )  Real 

SIGMA (25)  Real 


4.6 

.070 

.100 

.100 

3.0 

4.7 

.060 

.100 

.100 

3.0 

4.8 

.050 

.080 

.  005 

3.0 

4.9 

.040 

.070 

.005 

3.0 

5.0 

.035 

.055 

.002 

3.0 

0 

0 

0 

0 

0 

GLOSSARY  OF  PROGRAM  VARIABLES 


TOTMAS (25) 
TAME (25) 
ALTCHK ( 25 ) 
A ( 25 , 25) 

B ( 25 , 50 ) 

C  (25) 


Description 

Total  mass  or  activity  per  group 

Initial  lofted  altitude  for  each 

mono- size  group 

Mono-size  group  mean  radii 

Vertical  standard  deviation  for 

each  mono- size  group 

Total  mass  or  activity  per  sample 

Sample  collection  times 

Sample  altitudes 

Working  matrix 

Matrix  inversion  work  area 

Matrix  multiplication  solution 

vector 

3 

Fallout  density  (kg/m  ) 

2 

Acceleration  of  gravity  (m/s  ) 
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Variable 

Name 

Real/ Integer 

Description 

YKT 

Real 

Total  weapon  yield 

PLUTO 

Integer 

Number  of  samples 

MASS 

Real 

Total  mass  or  activity  lofted  by 

the  weapon 

I 

Integer 

Counter 

TONE 

Real 

Temporary  time  storage  location 

MARK 

Integer 

Counter 

BETA 

Real 

Distribution  log- slope  argument 

ALPHAl 

Real 

Starting  distribution  median 

radius 

ISTEP 

Real 

Median  radius  iteration 

increment 

ALPHA2 

Real 

Final  distribution  median  radius 

NUMGRP 

Real 

Number  of  mono-size  groups 

MOMENT 

Real 

Lognormal  distribution  moment 

RAT 

Real 

Temporary  storage  variable 

BP 

Real 

Distribution  log-slope 

PDINC 

Real 

Area  increment 

PDCHK 

Real 

Incremental  variable 

Z 

Real 

Working  variable 

PD 

Real 

Working  variable 

RIGHT 

Real 

Temporary  storage  variable 

LEFT 

Real 

Temporary  storage  variable 

FLAG 


Integer 


Counter 


Variable 

Name 


Real/Integer  Description 


COLT 

Integer 

Counter 

HOLT 

Real 

Temporary  storage  variable 

RAD 

Real 

Temporary  storage  variable 

INALT 

Real 

Temporary  storage  variable 

SIG 

Real 

Temporary  storage  variable 

IM 

Real 

Hopkins'  empirical  intercept 

SM 

Real 

Hopkins'  empirical  slope 

ID 

Real 

Conners'  empirical  intercept 

SD 

Real 

Conners'  empirical  slope 

DELTAZ 

Real 

Mono- size  group  vertical 

thickness  (m) 

LOGYKT 

Real 

Natural  log  of  YKT 

RHO 

Real 

Air  density 

AIDA 

Real 

Air  dynamic  viscosity 

LK 

Real 

Atmospheric  slope 

TZ 

Real 

Air  temperature 

PZ 

Real 

Air  pressure 

TEMPAL 

Real 

Temporary  storage  variable 

R2CD 

Real 

Reynolds  number  squared  times 

the  drag  coefficient 

REY 

Real 

Reynolds  number 

VZ 

Real 

Particle  fall  velocity  (m/s) 

SLPDRG 

Real 

Slip-drag  correction 

DELTAT 

Real 

Particle  fall  time  increment 

Variable 

Name 


Real/ Integer  Description 


+ 


J 

Integer 

T 

Integer 

• 

DUMMY 

Real 

ALTY 

Real 

SIGY 

Real 

# 

ALTZ 

Real 

DZ 

Real 

• 

SUMMER 

Real 

FIGS 

Integer 

DTHOLD 

Real 

• 

GAUSS 

Real 

M 

Integer 

• 

K 

Integer 

N 

Integer 

Q 

Integer 

• 

P 

Integer 

BF 

Real 

H 

Integer 

• 

F 

Integer 

U 

Integer 

TEMP 

Real 

TOTAL 

Real 

Sample  counter 
Group  counter 

Temporary  storage  variable 
Temporary  storage  variable 
Temporary  storage  variable 
Current  particle  altitude  (m) 
Distance  traveled  by  a  particle 
in  a  time  increment  DELTAT 
Time  of  fall  counter 
Flag 

Temporary  storage  variable 

Vertical  Gaussian  Value  at 

sample  altitude  at  sample  time 

Counter 

Counter 

Counter 

Counter 

Counter 

Temporary  storage  variable 

Counter 

Counter 

Counter 

Temporary  storage  variable 
Temporary  storage  variable 
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Variable 

Name 

X 

Y 

YOU 

WIN 

SEAL 

CRAB 

TAG 

MINUSE 

AVERAG 

FISH 

PIG 

OMEGA 

QUICK 

DENOM 

PARTY 

HOPKIN 

USAIR 


Real/ Integer 
Integer 
Integer 
Real 
Real 
Real 
Real 
Integer 
Integer 
Real 

Real 

Real 

Integer 

Integer 

Real 

Subroutine 


Subroutine 


Subroutine 


Description 

Counter 

Counter 

Temporary  storage  variable 
Temporary  storage  variable 
Temporary  storage  variable 
Temporary  storage  variable 
Flag 
Counter 

Average  of  solution  vector 
elements 

Standard  deviation  of  solution 
vector  elements 
Real  form  of  PLUTO 
Program  stop  flag 
New  SEARCH  parameter  flag 
Long/ short  search  flag 
Calculates  the  mean  radii  for 
the  equal  mass/activity  groups 
for  each  trial  distribution 
Calculates  the  stabilized 
altitudes  and  standard  deviations 
for  each  mono-size  group 
Calculates  RHO  and  AIDA  for  any 
given  altitude 
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Abstract 
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size  distribution  that  Is  positioned  In  the  atmosphere  by  empirical 
relationships  developed  by  Hopkins  and  Connors.  Davies -McDonald  fall 
meohanlos  are  used  to  model  the  falling  particles  In  the  nuolear  cloud. 

The  amount  of  mass  at  each  sample  altitude,  at  each  sample  time  Is 
calculated  from  the  cloud  model  and  compared  to  the  amount  of  mass  found 
in  the  actual  cloud  samples.  When  the  calculated  masses  equal  the  actual 
masses,  the  partiole  distribution  used  to  construct  the  stabilized  cloud 
is  the  correct  one.  A  computer  oode  for  this  numerical  analysis  is  also 
presented. 

The  computer  code  is  tested  using  hypothetical  filter  sample  data 
constructed  from  a  known  particle  size  distribution.  Additionally,  an 
input  parameter  sensitivity  analysis  Is  conducted. 

Actual  nuclear  cloud  sample  data  from  the  Redwing  series,  shot  ZUNI 
is  analyzed  using  this  numerical  method  of  airborne  nuclear  cloud  sample 
analysis.  The  outcome  of  the  ZUNI  sample  analysis  is  somewhat  inconclusive 
in  that  it  does  not  pinpoint  a  distribution.  However,  based  on  the  results 
of  the  model  sensitivity  analysis,  the  ZUNI  sample  analysis  Indicates  that 
the  particle  size  distribution  of  the  stabilized  ZUNI  cloud  may  be 
lognormal  with  a  log-slope  between  2.9  and  3»9  but  is  definitely  not  less 
than  2.7  nor  greater  than  5.0. 

Finally,  the  the  theory  for  an  alternative  method  of  airborne 
sample  analysis  is  presented.  Th3,s  method  uses  the  relative  number  of 
particles  of  each  size  found  in  airborne  samples to  unfold  the  stabilized 
nuclear  cloud  particle  size  distribution. 


